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Abstract 

We study Drell-Yan (DY) dilepton production in proton(deuterium)-nucleus and in 
nucleus-nucleus collisions within the light-cone color dipole formalism. This approach 
is especially suitable for predicting nuclear effects in the DY cross section for heavy ion 
collisions, as it provides the impact parameter dependence of nuclear shadowing and 
transverse momentum broadening, quantities that are not available from the standard 
parton model. For p{D) + A collisions we calculate nuclear shadowing and investigate 
nuclear modification of the DY transverse momentum distribution at RHIC and LHC 
for kinematics corresponding to coherence length much longer than the nuclear size. 
Calculations are performed separately for transversely and longitudinally polarized DY 
photons, and predictions are presented for the dilepton angular distribution. Further- 
more, we calculate nuclear broadening of the mean transverse momentum squared of 
DY dileptons as function of the nuclear mass number and energy. We also predict 
nuclear effects for the cross section of the DY process in heavy ion collisions. We found 
a substantial nuclear shadowing for valence quarks, stronger than for the sea. 

PACS: 13.85.Qk; 24.85. -hp; 24.70.-Fs 

Keywords: Drell-Yan process; nuclear shadowing; heavy-ion collisions 



1 



1 Introduction 



The cross section for the Drell-Yan (DY) process at the energies of the SPS suggests rather 
weak nuclear effects, if any at all (although measurements at small Feynman xp are the only 
data available). However, the fixed target experiment E772 at Fermilab at 800 GeV ^ shows 
a sizable nuclear suppression at large xp- Although this suppression results from a compli- 
cated interplay between energy loss and shadowing [0, ^, shadowing effects are expected to 
be much stronger and span the entire range oi xp at the energies of the Relativistic Heavy 
Ion Collider (RHIC) and the Large Hadron Collider (LHC). 

Relying on the standard parton model for proton-nucleus collisions, one can predict the 
DY cross section integrated over transverse momentum employing QCD factorization and 
data for the nuclear structure functions measured in deep-inelastic scattering (DIS)Q How- 
ever, no reliable means to calculate nuclear effects in the transverse momentum distribution 
within the parton model is available. Therefore, one has to parameterize nuclear broadening 
according to Brownian motion in the plane of transverse momentum and fit to data (see 
for a review). Perhaps, this situation will be improved by higher-twist factorization theorems 
0] in the future. 

Moreover, parton model predictions are doubtful even for the integrated DY cross section 
in nucleus-nucleus collisions. Indeed, compared to pA collisions, this case requires knowledge 
of the impact parameter dependence of nuclear shadowing and, of course, the introduction 
of an additional integration over the impact parameter [see Eq. (|^)]. Neither DIS nor the 
DY reaction on nuclei provides the information needed for this. In principle one can access 
such information relying, for instance, on knowledge of the number of so-called gray tracks 
and using simple cascade models. However, this possibility has never been realized either for 
the DIS or the DY process. In view of this problem, it might be tempting to assume, as in a 
recent analysis |p , that nuclear shadowing is independent of impact parameter. Clearly, this 
cannot be correct and can lead to unphysical results such as a DY cross section in heavy ion 
collisions independent of centrality. It is known, however, for many processes, that peripheral 
collisions are similar to the free A^A^ interaction, while central collisions should manifest the 
strongest nuclear effects. 

In this paper we calculate nuclear shadowing for the DY cross section using the light- 
cone (LC) dipole approach suggested in 0. This approach formulates the DY reaction in 
the rest frame of the target, which allows one to overcome these problems in a simple way by 
taking advantage of the description of the DY process as bremsstrahlung of a heavy photon 
from a beam quark that subsequently decays into the lepton pair as in fig. |1]. Although 
this looks very different from the more familiar DY mechanism which interprets the 



DY reaction as quark-antiquark annihilation, it is known that the space-time interpretation 
of high-energy reactions is not Lorentz invariant and depends on the frame of reference. 
Indeed, it is only in a fast-moving frame where Feynman's picture of the colliding particles 



^Thc analysis of data Q based on the DGLAP evolution equations still neglects effects of saturation Q 
that should be important once shadowing sets in. Additionally, no data for DIS on nuclei are available for 
small Bjorken x relevant for LHC. Existing data for DIS and DY are not sensitive to the amount of gluon 
shadowing, which had to be parameterized ad hoc in 0]. 
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as bundles of non-interacting partons with no (or small) transverse momenta is valid and 
the DY reaction can be formulated consistently in terms of the parton density of the proton 
(or nucleus). 




Figure 1: In the target rest frame, DY dilepton production looks like 
hremsstrahlung . A quark or an anti-quark from the projectile hadron 
scatters off the target color field (denoted by the shaded circles) and ra- 
diates a massive photon, which subsequently decays into the lepton pair. 
The photon decay is not shown. The photon can be radiated before or 
after the quark scatters. 

In the target rest frame, a quark of the incident hadron fluctuates into a state that 
contains a massive photon (dilepton) and a quark. Interaction with the target breaks down 
the coherence of the fluctuation and the 7* is freed. Correspondingly, the cross section of 
the process — >■ 7* X has a factorized form ||^, |Tl], |12|, |13|, |I4| , 



din a 



where \l/^.g(a,p) is the LC distribution amplitude in Eqs. ( |B.2| ) or ( |B.3| ) for having a quark- 
photon (transversely or longitudinally polarized) fluctuation with transverse separation p 
and relative fractions a and 1 — a of light-cone momenta carried by the photon and quark, 
respectively. For the DY reaction in pp-scattering, the dipole cross section needed in (|l|) is 
the same as in DIS off a proton Q Note, that in the two graphs for bremsstrahlung, fig. ^ 
the quark scatters at different impact parameters, depending on whether it scatters when in 
the |7*g)-state (right) or not (left). This leads to the appearance of the dipole cross section 
agq{ap,X2) in (|l|), although there is actually no physical dipole in this process The 



light-cone Fock-state formulation of QCD is discussed in detail in the recent review |[T5|] . 
We use standard kinematical variables, 

Xi = , X2 = , (2) 

s s 

with Xi — X2 = xp and X1X2 = (M^ -|- g|,) / s, where Pi, P2 and q are the four-momenta of the 
beam, target and the virtual photon, respectively; = and are the dilepton invariant 
mass squared and transverse momentum, respectively; s = {Pi + P2)^- 

The frame dependence of the space-time interpretation of the DY process can be illus- 
trated by the different meanings of xi in different reference frames: it is well known that in 

^This can be proven in leading-log order or otherwise justified by referring to QCD factorization. 
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the standard picture of DY ||T^, Xi is the momentum fraction of the projectile quark anni- 



hilating with the target antiquark. However, evaluating the scalar product (Q) in the target 
rest frame shows that the projectile quark carries momentum fraction x = Xi/a > Xi of its 
parent hadron and, correspondingly, Xi is the momentum fraction of the proton taken away 
by the photon. This is not a contradiction, since the projectile quarks in the two reference 
frames are different particles. 

In the case of a pA-scattering, one has to distinguish between two limiting kinematical 
regimes for the DY reaction. On the one hand, there is the regime of short coherence time tc, 
which can be interpreted as the mean fluctuation lifetime. If tc is much shorter than the mean 
internucleon separation no effect of coherence (shadowing) is expected. On the other hand, 
in the regime of long coherence time compared to the nuclear radius, tc ^ Ra, interference 
of the multiple interaction amplitudes with bound nucleons will affect the probability of 
breaking down the coherence of the fluctuation and releasing the dilepton on mass shell, 
i.e. shadowing (sometimes antishadowing) occurs. These interferences are controlled by the 
longitudinal momentum transfer Qc = 1/tcin the process of 7* radiation by a projectile quark 
of energy Eg, qN ^ qX, 

M^'n - ml 

where we assume energy conservation, and the invariant mass squared of the 7*g pair is 

M"^ "^a Qt 



One arrives at a similar estimate with help of the uncertainty principle. Indeed, a quark 
can violate energy conservation by fluctuating into 7*g for a time At ~ l/(A£y*q — niq) 
in the quark's rest frame. Applying the Lorentz gamma-factor 7 ~ '2Eq/ {M^*q + m^) one 
reproduces the lifetime tc = 1/qc 'v^ the lab frame as given by (H). 

The intuitive space-time pattern related to the coherence time for DY pair production off 
nuclei is rather obvious. In the limit of short coherence time (relevant for the SPS energy) 
the initial state interactions are predominantly soft, since the hard fluctuation containing 
the heavy dilepton appears only deep inside the nucleus and is immediately freed on mass 
shell via the interaction with a bound nucleon. 

On the other hand, if the coherence length substantially exceeds the size of the nucleus (as 
expected for the energies of RHIC for xp greater than about 0.5 and the LHC for essentially 
all xf)i the hard fluctuation is created long in advance of the interaction with the nucleus, 
which acts as a whole in freeing the fluctuation. Since different target nucleons compete 
with each other, the DY cross section is subject to shadowing. 

In this paper, we study nuclear effects in the limit of long coherence time. This is a most 
interesting regime, where interference effects are maximal and all the nucleons having the 
same impact parameter participate coherently in the DY process. A special advantage of 
the color-dipole approach is that it allows one to incorporate nuclear shadowing via a simple 
eikonalization of the dipole cross section HTBl 0] (see next section). This follows from the 



fact that in this limit the dipole size is "frozen" by Lorentz time dilation. 
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A projectile quark can develop more complicated fluctuations that involve gluons in addi- 
tion to the heavy photon. These correspond to Fock states |g7*(j>'), 157*2(7), etc. Interaction 
of these fluctuations with the nucleus is also affected by shadowing, which may be even 
stronger than for \q'j*) provided that the fluctuation lifetime is long compared with nuclear 
size. This additional shadowing is, in the parton model, related to the shadowing of gluons 
resulting from gluon fusion at small X2 (see Sect. ||). 

One might think that in the case of pA collisions shadowing for DY can be easily predicted 
relying on QCD factorization and using data for shadowing in DIS off nuclei. However, data 
at small x are available only at very low Q^, where neither factorization nor DGLAP evolution 
are expected to be valid. Additionally, one should be cautious applying factorization at 
large xi where, as we pointed out above, higher twist corrections are rather large [^. In 
particular, the Bjorken x of the target, X2, reaches its minimal value as xi — *■ 1; therefore, 
factorization predicts the maximal strength for shadowing. However, shadowing for the DY 
process vanishes in this limit. Indeed, since a > xi, the invariant mass Eq.(^) increases (for 
massive quarks and/or nonzero qt), leading to the disappearance of the coherence length. 

We do not make a fit to the observed shadowing in DIS on nuclei. We rather follow 
the logic of the conventional Glauber approach, where one is permitted to make fits to the 
data for nucleon-nucleon collisions but then predict nuclear effects in a parameter-free way. 



Indeed, we use the phenomenological dipole cross section on a nucleon target [T^ fitted 
to the data for ep DIS from HERA that cover a range of energies much higher than that 
available from fixed-target nuclear shadowing data. The DY cross section calculated with 
this phenomenological cross section is supposed to include all higher-order corrections and 
higher twist effects. 

Another important advantage of our approach is the possibility of calculating nuclear 
effects in the transverse momentum distribution of DY pairs. This is a difficult problem 
within the parton model. The phenomenon of nuclear broadening of the dilepton transverse 
momentum looks very different at low (short tc) and high (long tc) energies. If tc is short, a 
hard fluctuation containing the heavy dilepton is created deep inside the nucleus just before 
the interaction that releases it. Meanwhile, the incident hadron may have soft initial state 
interactions in the nucleus. Although these do not generate DY dileptons, they do increase 
the mean transverse momentum of the fast partons of the projectile. Indeed, a fast parton 
experiencing multiple interactions performs something like Brownian motion in the plane of 
transverse momenta. As a result, the parton arrives at the point of the DY pair creation with 
increased transverse momentum. Information about the enhanced transverse momentum of 
the projectile quark in a nucleus compared to a proton target is carried away undisturbed by 
the dilepton pair. Nuclear broadening in the limit of short coherence time was investigated 
previously in 



At first glance, the observed broadening of the dilepton transverse momentum in the 
parton model might be viewed as arising from an increased transverse momentum of quarks 
and antiquarks in the nucleus. However, such a conclusion contradicts the usual picture of 
a nucleus boosted to the infinite momentum frame. There, nucleons and their parton clouds 
are well separated and do not overlap at large X2, just as in the nuclear rest frame. We 
know that only at very small X2 do the parton clouds overlap and fuse leading to nuclear 



5 



shadowing. Such a fusion process results not only in suppressing the parton density, but it 
also increases the transverse momenta of the partons. Thus, shadowing and qx broadening at 
small X2 are closely related processes; no broadening is possible without shadowing. However, 
the regime of short tc corresponds to large X2, where neither shadowing, nor qx broadening 
is expected for the nuclear parton distribution function. Thus, we face a puzzle: nuclear 
broadening of the transverse momentum distribution of DY pairs calculated in the nuclear 
rest frame and observed experimentally has no analog within the parton model. This puzzle 
was resolved a long time ago by Bodwin, Brodsky and Lepage [l^, who found that in the 



regime of short coherence length initial state interactions leading to qx broadening violate 
QCD factorization and should not be translated into a modification of the nuclear quark 
distribution function. Likewise, initial state energy loss ||^, ^ cannot be translated to a 
modification of x-distribution of partons in nuclei. This is a general statement which is 
applied to other hard reactions like high-g^ hadron production, etc. 

In the regime of long coherence time tc ^ Ra relevant to RHIC and LHC, a very different 
mechanism is responsible for broadening of the transverse momentum distribution |]12|. A 
high-energy projectile quark emits a dilepton fluctuation (via a virtual time-like photon) long 
before its interaction with the nucleus. These components of the fluctuation, the recoil quark 
and the dilepton, do not "talk" to each other because of Lorentz time dilation. Therefore, 
multiple interactions of the quark in nuclear matter seem to have no further influence on the 
produced dilepton and no broadening of the transverse momentum is expected. However, 
this conclusion is not correct. While it is true that different ingredients of the fluctuation 
cannot communicate, not all fluctuations contribute to DY pair production in the same way: 
The harder the fluctuation, that is the larger the intrinsic relative transverse momentum 
between the quark and dilepton, the stronger the kick from the target required for loss of 
coherence, i.e. for the fluctuation to be produced on mass shell. Large 7*g-fluctuations 
(which have a low intrinsic transverse momentum) have a high probability to lose coherence 
in the first scattering on the surface of the nucleus. Therefore, the contribution to the DY 
cross section from these configurations scales as A^/^, resulting in a nuclear suppression of 
the DY cross section at low qx- Due to multiple interaction of the projectile quark a nucleus 
provides a larger transverse momentum than a nucleon target and hence is able to break up 
smaller fluctuations than a nucleon. This enhances the DY cross section at intermediate qx- 
At high qx multiple scatterings are not important any more, just a single high qx scattering 
dominates, eliminating the nuclear enhancement. This is a result of the power qx dependence 
of the cross section predicted by QCD. Since small qx dileptons are suppressed while large 
qx dileptons are not, the mean transverse momentum of DY pairs is increased. Thus, at long 
tc, the mechanism behind nuclear broadening is color filtering pO| . 

This paper is organized as follows. We explain the main ideas of the light-cone approach in 
the introduction. The key ingredient of this method, the universal color dipole cross section 
for a qq pair interacting with a nucleon, is known from phenomenology. In Sect. |^, we explain 
how nuclear effects are treated in the color dipole approach. In particular, we describe, how 
nuclear gluon shadowing has to be included along with the gg-nucleus cross section. The 
results for gluon shadowing are presented in Sect. ^. The results of our calculations for 
the DY cross section in p{D) + A collisions, and predictions for RHIC and LHC, can be 
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found in Sect. ^. Nuclear modification of the DY pair transverse momentum is calculated 
for the energies of RHIC and LHC in Sect. ^ We find that the so called Cronin effect, 
nuclear enhancement of the DY cross section at medium-large g^, is nearly compensated 
at RHIC energies but is expected to have a large magnitude at LHC. We conclude that 
nuclear broadening of the DY transverse momentum squared diverges logarithmically for 
transversely polarized photons if nuclear shadowing occurs. Differences in nuclear effects 
for radiation of longitudinally and transversely polarized photons lead to a specific nuclear 
modification of the DY polarization. Corresponding predictions are presented in Sect. |[ In 
addition, we find substantial deviation from the so called Lam- Tung relation ||21|]. Indeed, 
this relation is not supported by data, which is difficult to explain within the standard parton 
approach. In Sect. |^ we address the more difficult problem of nuclear effects in heavy ion 
collisions. We follow conventional wisdom and simplify the problem by employing QCD 
factorization. Nuclear shadowing for both sea and valence quarks is calculated within the 
LC dipole approach. In contrast to usual expectations, shadowing for valence quarks turns 
out to be larger than for the sea. We summarize the results and observations of this paper 
and present an outlook for further development and application of the LC dipole approach 
in Sect. |8|. 



2 The gg-nucleus cross section 



In order to calculate the DY cross section in pA scattering, one has to replace a^^ in Eq. ([1|) 
by the color dipole cross section on a nucleus, (jfq{p, x), which is easy to calculate within the 
color dipole approach. 

In the limit of long coherence time, the projectile quark may be decomposed into a series 
of Fock-states with frozen transverse separations. Since partonic configurations with fixed 
transverse separations in impact parameter space are interaction eigenstates [0, cr^(p, x) 



may be calculated using Glauber theory ||2^, i.e. via simple eikonalization of the gg-nucleon 
cross section, 

A' 

qq\ 



{p, 



1 



— a>,a;)T^(6) 



(5) 



where h is the impact parameter, A is the nuclear mass number and 

TA{b) = / dzpA{b, z) 



(6) 



is the nuclear thickness, i.e. the integral over the nuclear density. The possibility of eikon- 
alization of a color dipole propagating through a nucleus suggested in [0 has been more 
rigorously justified recently in [jl8|, |2^, |25|. We mark o"^ in Eq. (1) with a tilde to indicate 
that it misses important contributions from higher Fock components, which motivates us to 
introduce an improved eikonal formula, Eq. ([T3|) , below. Note that at larger X2 ^ 0.01, when 
the coherence length becomes shorter than the nuclear radius, the eikonal approximation 



breaks down and one has to employ the Green function technique developed in [12 . 
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Figure 2: The eikonal formula (Qj takes only multiple rescatterings of the 
\qq)-Fock component into account. This figure illustrates the amplitude 
for double scattering (left). When the amplitude is squared (middle), the 
gluon rungs combine into gluon ladders (Pomerons), which are enclosed 
into each other. In Regge theory, this contribution to the cross section is 
expressed in terms of the Pomeron-Pomeron-Reggeon vertex (right). 



The single-scattering term can be obtained by expanding (|]) to first order in agg{p,x). 
The dipole interacts with the target by exchange of a gluonic colorless system, the so called 
Pomeron. The unitarity cut of such an amplitude reveals multiple gluon radiation that is 
related to higher Fock states within the LC approach in the target rest frame. Thus, for 
single scattering, cr^ takes all Fock states of the projectile parton into account, not only 
\qq). The energy dependence of the dipole cross section is generated by the phase space of 
gluons from higher Fock states \qqG), \qqGG), . . .. Indeed, in the Born approximation, i.e. 
two gluon exchange, would be independent of x. 

Calculation of cr^ from first principles is still a challenge. We rely on phenomenology and 
employ the parameterization of Golec-Biernat and Wiisthoff Jl^ motivated by the saturation 
model. 



exp 



A{x/xo)^ J 



(7) 



where Qo = 1 GeV and the three fitted parameters are ao = 23.03 mb, xq = 0.0003, and 
A = 0.288. This dipole cross section vanishes oc at small distances, as implied by color 
transparency, and levels off exponentially at large separations. The authors of |T^ are able 
to fit all available HERA data with a quite low and can also describe diffractive HERA 
data. 

We now turn to the multiple scattering terms. Describing shadowing for DY is simplified 
if we make use of the fact that the dipole cross section entering the formula for dilepton 
production, Eq. (p, is the same as that needed to calculate the DIS cross section. We may 
thus illustrate the physics of Eq. @ in fig. ^, where the double-scattering term for a qq- 
dipole is depicted. In terms of Regge phenomenology, the double scattering of the gg-pair 
corresponds to the Pomeron-Pomeron-Reggeon vertex. Note that (|^) accounts for not only 
the double-scattering term, but also for all higher-order rescatterings of the gg-pair. The 
ra-fold scattering graph has n gluon ladders, which are enclosed into each other. 

Rescatterings of higher Fock states, containing gluons are omitted in (|^). This equation 
corresponds to gluon radiation in the Bethe-Heitler regime, where each interaction of the 
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Figure 3: At high energy, the lifetime of higher Fock states becomes long 
enough for multiple scattering. Shown here is the double-scattering am- 
plitude for the \qqG)-Fock state. In Regge theory, this process is expressed 
in terms of the triple- Pomeron vertex (right). The eikonal formula ([^ is 
improved by multiplying aqq by the gluon shadowing ratio Rq, Eq. 
to include these contributions also. 



quarks leads to independent radiation of the entire spectrum of gluons. However, interfer- 
ences known as the Landau-Pomeranchuk effect, lead to a suppression of gluon radiation. 
This shadowing correction is taken into account via multiple interactions of radiated gluons 



in the nucleus and leads to an additional suppression of the DY cross section. At low 
X the lifetime of these higher Fock states becomes significantly longer than the mean inter- 
nucleon distance, and they can scatter more than once inside the nucleus, as illustrated in 
fig. 1^. In this case, which occurs at RHIC and LHC energies, (H) needs to be modified to 
include also these rescattering. 

In order to include processes like the one illustrated in fig. |^, it is useful to note that 
the rescattered gluon can be interpreted as the first rung of a single gluon ladder exchanged 
between the gg-pair and the target. In Regge phenomenology, rescattering of gluons leads 
to the triple- Pomeron vertex fig. ^( right), which can be regarded as a correction to the 
single-scattering term. More precisely, it leads to a reduction of the nuclear gluon density 
Ga because the two Pomerons from the target in fig. ^(right) fuse into a single one before 
interacting with the pair. Indeed, multiple scatterings of higher Fock states containing 
gluons are known as the effect of gluon shadowing and lead to an additional suppression 



of the DY cross section. In the infinite momentum frame of the nucleus, gluon clouds of 
different nucleons overlap and fuse at small x ^, fig. ^( right), thereby reducing the gluon 



density at small x. Although the corrections for gluon rescatterings the nuclear rest frame 
and gluon fusion (in the nuclear infinite momentum frame) look very different, this is the 
same phenomenon seen from different reference frames. Of course, observables are Lorentz 
invariant, and both effects lead to a reduction of the DY cross section. 

Thus, the Pomeron-Pomeron fusion process in fig. |^(right) can be taken into account by 
multiplying o"^(p, x) by the gluon shadowing ratio 

Rc{x, Q') = ^^^^^ - 1 - ARaix, Q% (8) 
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as calculated in Appendix A. The single-scattering term then reads 



<(x,p) = a^'Jp^x) I d;'bRGix,Q-',b)TA{b) + ..., 



1-j I d%ARG{x,Q\b)TAib) 



+ ... 



(9) 



where gluon shadowing as a function of impact parameter b is given in Eqs. ( A.2CI|) and 



The first term in the brackets of stands for the direct exchange of a Pomeron, while the 
second (negative) term represents the correction due to the Pomeron fusion process depicted 
in fig. H This recipe becomes even more clear from the relation, valid at p — 0, El 



■qq 



p->0 



(10) 



In (|^), the gluon density of the proton in (|T0|) was replaced by the average gluon density of 
the nucleus, 



TT 



-a. 



.2 iP j'^^[''^-2 



11) 



A 

.p^y ■ ^ VP 

It is clear from (^ that the effective dipole cross section on a bound nucleon appears to be 
reduced due to gluon shadowing. We also see from (p^Of) that Rq has to be evaluated at a 
scale Q"^ = X/p^. 

The specific fusion processes that are included in (H) depends on the approximation in 
which gluon shadowing Rq is evaluated. For our actual calculations (see Sect. ^ nuclear 
shadowing for gluons is calculated within the Green function formalism for a \qqG) fluc- 



tuation propagating through nuclear medium as developed in |^9[. This means, that the 
single-scattering term in (^) is corrected not only for the 2IP IP Pomeron fusion term 
depicted in fig. |^(right), but also for all the ralP TP fusion processes. Moreover, the Green 
function approach properly describes the finite lifetime of the |ggG)-state. This is important, 
because even when the gg-fiuct nation lives much longer than the nuclear radius, the lifetime 
of the |gg(j')-state will be shorter. 

For the rescattering terms, we can account for higher Fock states in the same way as in 



namely by the replacement [pOl pll, p2] 



a^^{p,x)^a^^{p,x)RG{x,Q\b) , 
i.e. the improved formula Eq. (^) for the gg-nucleus section reads. 



afJp,x)=2 / d'b 



1 - ( 1 - ^ <(P, x)Rg{x, \/p\ b)TA{b) 



(13) 



This expression includes also the contribution of higher Fock states containing more than 
one gluon. The higher-order multiple interactions of the \qqG) Fock state correspond, as 
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Figure 4: The figure on the left shows the qq-nucleus cross section j^TSj ) 
divided by the nuclear mass number A for two different values of x. The 
two lower curves (solid and dashed) are calculated for gold (A = 197^ 
and the two upper curves for a proton. The figure on the right shows the 
qq-nucleus cross section divided by A times the dipole cross section (Qj. 
While large separations are strongly suppressed, small size dipoles are 
much less affected by the nucleus. Nuclear gluon shadowing is included 
in the calculation, as explained in the text. It slowly vanishes at small qq 
separations, which correspond to high . 



was mentioned, to multi-Pomeron fusion, nJP IP, while the Reggeon diagrams with 
nW mW (m > 2) are missing. Those diagrams should be incorporated via the Fock 
components \qqmG) containing two or more gluons. The modified expression Eq. (0) sums 
multiple interactions of the qq pair via mJP exchange (summed over m), each of which 
has the form of a fan nIP — s> IP (summed over n). We assume that each gluon in the 
Fock state \qqmG) experiences multiple interactions independently of other gluons. This 
assumption corresponds to the Gribov's interpretation |Q of Glauber eikonal shadowing; 
namely, the unitarity cut of an n-fold scattering term must contain a simultaneous cut of 
all n Pomerons. It therefore corresponds to a nonplanar graph describing the independent 
multiple interactions of n projectile partons (see discussion in [PT|). 

Note that the GLR equation ||^ can be obtained from Eq. (|13j) by expanding the r.h.s. 
up to quadratic order in a^. Though the higher order terms in this expansion may not be 
found in the large N^, leading log(l/x) limit of QCD their numerical infiuence on the 
DY cross section will be small because of the large dilepton mass. 

We can now proceed to calculate (T^(x, p) according to (|13]). We briefiy summarize our 
calculation of gluon shadowing {Rg) in section ^. The results for the nuclear dipole cross 
section cr^(p, x) are depicted in fig. ^. Since (|T^ is a high-energy approximation, valid when 
the lifetime of the gg-pair exceeds the nuclear radius, these results are relevant at RHIC 
and LHC energies. At lower energies however, one has to take transitions between different 
eigenstates into account . 



The plot on the l.h.s. of fig. |J shows the dipole cross section itself. First, we discuss the 
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gg-proton cross section. The two upper curves show this quantity for two different values of 
X typical for RHIC and LHC After a quadratic rise, o"^ levels off and assumes an energy- 
independent saturation value of 23.03 mb. The onset of saturation, i.e. the flattening of the 
dipole cross section as function of p, is controlled by the saturation radius. 



R'(^) = ^ -) > (14) 



Ql \x 

which decreases with energy. The energy dependence of correlates with p. At p -C i?^, 
the dipole cross section grows with energy with the hard Pomeron intercept A = 0.288, while 
at large separations, p ^ Rg, (0) becomes independent of energy. For more discussion on 
(0), refer to the original work [|17 



We now turn our attention to the qq nucleus cross section (T^(p, x), Eq. (|T3|), which 
is shown by the two lower curves in fig. ^ (left). In all calculations, we employ realistic 
parameterizations of nuclear densities from |^. In addition to the expected suppression 
due to nuclear shadowing, one also sees that the saturation value of cr^, which is almost at 
its geometrical limit 27iR\, is energy dependent. Moreover, (J^(p oo, x) is a decreasing 



function of energy. This is a consequence of the gluon shadowing in (p!3|). At very small x 



gluon shadowing becomes strong, (see fig. ^ and the qq nucleus cross section lies below its 
geometrical limit. The stronger the gluon shadowing, the smaller the saturation value of o"^. 

However, Rq — > 1 at p — > since ~ 1/p^ — oo (see fig. ||). 

The nuclear suppression of the dipole cross section is plotted on the right-hand side of 
fig. ^. Small sizes are less affected by the nucleus than large sizes. This illustrates the 
effect of color filtering ||20| , which is the mechanism behind nuclear broadening of transverse 
momenta (see section^. While small gg-pairs, which have large intrinsic transverse momenta 
according to the uncertainty principle, propagate through the nucleus almost undisturbed, 
large pairs (small transverse momenta) are absorbed, i.e. the coherence of the fluctuation 
is disturbed and the 7* is freed. This effect leads to an increase of the mean transverse 
momentum. 



3 Gluon shadowing in nuclei 

The nuclear shadowing for gluons needed as input for the gg-nucleus cross section in (IT 



is calculated in the LC Green function approach developed in [29|, where gluon shadowing 



is calculated from shadowing of the IggG) Fock component of a longitudinally polarized 



photon. In this section, we briefly review the approach of p9[ and present the results of our 
calculation for gluon shadowing as function of x, and the length L of the path in the 
nuclear medium. 

Longitudinal photons can serve to measure the gluon density because they effectively cou- 
ple to color-octet-octet dipoles. This can be understood in the following way: the light-cone 
wave function for the transition 72 — > gg does not allow for large, aligned jet configurations. 
Thus, unlike the transverse case, all gg dipoles from longitudinal photons have size 
and the double-scattering term vanishes like oc The leading-twist contribution for the 
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Figure 5: Gluon shadowing vs. the length of the nuclear medium L = 
2-^/^^"^^^, where b is the impact parameter and Ra the nuclear radius. 
All curves are for = 20 GeV^ but for different values of x. 



shadowing of longitudinal photons arises from the \qqG) Fock state of the photon. Here 
again, the distance between the q and the q is of order l/Q^, but the gluon can propagate 
relatively far from the gg-pair. In addition, after radiation of the gluon, the pair is in an 
octet state. Therefore, the entire |ggG)-system appears as a GG-dipole, and the shadowing 
correction to the longitudinal cross section is just the gluon shadowing we want to calculate. 

A critical issue for determining the magnitude of gluon shadowing is the distance the 
gluon can propagate from the gg-pair in impact parameter space, i.e. knowing how large the 
GG dipole can become. This value could be extracted from single diffraction data in hadronic 
collisions in 123] because these data allow the diffractive gluon radiation (the triple- Pomeron 
contribution in Regge phenomenology) to be unambiguously singled out. The diffraction 
cross section (oc p^) is even more sensitive to the dipole size than the total cross section 
(oc p^) and is therefore a sensitive probe of the mean transverse separation. It was found 
in ||29| that the mean dipole size must be of the order of tq = 0.3 fm, considerably smaller 
than a light hadron. A rather small gluon cloud of this size surrounding the valence quarks 
is the only way that is known to resolve the long-standing problem of the small small size 
of the triple- Pomeron coupling. The smallness of the GG dipole is incorporated into the LC 
approach by a nonperturbative interaction between the gluons. 
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Figure 6: The x- and Q"^ -dependence of gluon shadowing for carbon, 
copper and gold. The x-dependence is shown for = 20 GeV^ , while 
the figure on the right is calculated for x = 10~^. 



Note that the small value of tq dictated by data for diffraction is consistent with the 
results of other approaches. Indeed, the same small size characterizing gluonic fluctuations 



was found in the instanton liquid model in the QCD sum rule analysis of the gluonic 



formfactor of the proton and it also follows from lattice calculations |3^. Note that the 
value of To also limits the Q^-range where the approximation qqG ~ GG is valid. One has 
to ensure that ^ V'^O' otherwise the qq pair is not pointlike compared to the size of the 
entire Fock state. 

Our results for gluon shadowing as a function of the length of the nuclear medium at 
impact parameter b are shown in fig. |^. The calculations are performed for lead with a 
uniform nuclear density of pa = 0.16 fm~'^. Details are presented in [Appendix A| . The small 
size of the GG dipole leads to a rather weak gluon shadowing (except for specific reactions 
where the qq pair is colorless [0). For most values of x, gluon shadowing increases as a 
function of L as one would expect. At the largest value of x = 0.01, however, gluon shadowing 
becomes smaller as L increases, and Rq approaches 1. Although this behavior seems to be 
counterintuitive, it can be easily understood by noting that at x = 0.01 the coherence length 
of the \qqG)-Fock state becomes very small and the formfactor of the nucleus suppresses 
shadowing [RU| . The curves shown in fig. H are the ones which actually enter our calculation 



for DY via the gg-nucleus cross section (13). The values of x entering our calculation are 
X ^ 10-3 for RHIC and x ^ 10"^ for LHC. 

We also calculate gluon shadowing as function of x at fixed and as a function of 
at fixed x, integrated over the impact parameter b. The results are shown in fig. ^. In 
the left-hand plot, one observes that gluon shadowing vanishes for x > 0.01. This happens 
because the lifetime of the IggG) -fluctuation becomes smaller than the mean internucleon 
distance of ~ 2 fm as x exceeds 0.01. Indeed, in [RUI an average coherence length of slightly 
less than 2 fm was found for the \qqG)-state at x = 0.01 and large ^ l/''"o- Note that 
gluon shadowing sets in at a smaller value of x than quark shadowing because the mass of 
a |ggG')-state is larger than the mass of a |gg)-state. This delayed onset of gluon shadowing 
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was already found in |2^. We also point out that gluon shadowing is even weaker than quark 
shadowing in the x-range plotted, because the small size of the GG-dipole overcompensates 
the Casimir factor in the GG-proton cross section, a^g. = (9/4) cr^. The plot on the right- 
hand side of fig. ^ shows the Q^-dependence of gluon shadowing and clearly demonstrates 
that gluon shadowing is a leading-twist effect, with Rg only very slowly (logarithmically) 
approaching unity as — * oo. 



4 Nuclear shadowing for DY pair production in pA- 
and iDA-collisions 

Nuclear shadowing for the DY process was first observed in proton-nucleus {pA) collisions at 
large Xi? in the E772 experiment. The shadowing effect will also be present in the energy 
range of RHIC and LHC Since RHIC will probably first measure the DY cross section from 
deuterium-nucleus {DA) rather than pA collisions, we perform calculations for both, pA and 
DA collisions. 

The dipole formulation provides the following explanation of shadowing in the DY pro- 
cess. When the coherence time is long, one of the projectile quarks develops a 7*g-fluctuation 
long before it reaches the target. If the transverse momentum transfer from the target is 
large enough to resolve the fluctuation, the virtual photon is freed and eventually is observed 
as a lepton pair in the detector. In the case of a nuclear target, the set of struck nucleons 
compete to free the virtual photon. If the |7*g)-state has a very small transverse size, it can 
propagate through the entire nucleus because none of the bound nucleons can provide a kick 
strong enough to resolve the |7*g) structure in the incident quark. These small fluctuations 
have the same small probability to interact with any nucleon, so they will not be shadowed. 
On the other hand, if the fluctuation is sufficiently large in size, only a small momentum 
transfer is necessary to resolve the photon. Thus, the coherence of a large fluctuation will 
be destroyed with high probability already in the first collision on the surface of the nucleus. 
Nucleons deeper in the nucleus do not add much to the probability of freeing the 7*. Thus, 
this probability nearly saturates for these extremely large fluctuations, and the DY cross 
section will scale like A^/^. From these considerations, we can find two necessary conditions 
for shadowing PP|, 



The 7*g fluctuation must have a lifetime long enough to allow for at least two scatter- 
ings during the coherence time tc. 

The 7*g fluctuation must have a large freeing cross section^, i.e. its transverse size must 
be sufficiently large. 



•^One should distinguish between the freeing and the total cross sections of a fluctuation. The latter 
is always large for a colored quark and all its fluctuations, while the former is driven by the difference 
between scattering amplitudes of different fluctuations and |(77*) in our case). Since it is only the quark 
that interacts in each of these Fock states, the freeing cross section is controlled by the relative transverse 
displacement of the quarks within different fluctuations. This is how the dipole cross section comes about. 
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Figure 7: Shadowing for the total DY cross section in proton - gold (upper 
curves) and deuterium - gold (lower curves) collisions at the energies 
of RHIC and LHC as function of Feynman xp and dilepton mass M"^, 
respectively. The left-hand figure is calculated for M = 4.5 GeV. The 
figure on the right for xp = 0.5. 



The first condition is assumed to be fulfilled throughout this paper, where we consider only 
the case of infinite tc- The dependence of shadowing on the fluctuation size is encoded in 
the gg-nucleus cross section (|13|). 

Note that since the 7*g-fluctuation is formed long before it reaches the the target, the 
dilepton is unaffected by quark energy loss. Thus, the entire suppression of the DY cross 
section at very low X2, say X2 < 0.001, is due to shadowing, and we do not need to worry 
about energy loss. This is different at the lower fixed-target energies |]1]], where the observed 
depletion of the DY cross section originates from a combination of shadowing and energy 
loss [0, ^. The complimentary behavior of shadowing and energy loss is discussed in more 
detail in : long tc means that only shadowing occurs, while for short tc, one observes only 
energy loss. 

To obtain the cross section for an incident hadron, the partonic cross section (|l]) has to 
be weighted with the quark (and antiquark) distributions, q{x), of the projectile hadron, and 
one has to include in addition the factors necessary to account for the decay of the virtual 
photon into the dilepton. For an incident proton, this cross section becomes 

d^a{pA — > IIX) OLem ^\ da -sr-^ 2 f /^l A4'2\ /^i 



dM^dxp 37rM2xi + X2 



d(y{qA 7*X) 



d\\ia ' 

where Zf is the charge of a quark of flavor /. We assume that the the same expression ([T5| ) 
applies for both proton and deuteron projectiles, so that the only difference between these 
cases is that the flavor sum ranges over the quarks of the proton and neutron in the case 
of an incident deuteron. Nuclear effects in the deuterium structure function and finite-size 
effects are neglected, and isospin symmetry is assumed. 
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For a calculation that can actually be compared to data, we employ the CTEQ5L pa- 
rameterization (taken from CERNLIB |^) for q^j. Note that the projectile quark 
distributions enter at large x = Xi/a > xi, where they are well known. Thus, the uncer- 
tainty arising from the choice of parton distributions is minimal. However, these parton 
distributions are different for the proton and deuterium, so that the pA and DA DY cross 
section are not trivially related. Shadowing can now be obtained by evaluating (^) and (|l]), 
with cr^(p, X2) taken from ([131) , and dividing by A times the analogous calculation with the 
gg-proton cross section (|^). In the case of the deuterium projectile, we divide by 2A. The 
nuclear density parameterizations are taken from ||35[| . 

The result as function of xp and dilepton mass M at different energies is shown in fig. ^ 
For each energy (RHIC and LHC) we calculated pAu and DAu collisions and normalized 
both to pp collisions. Note that the DAu curve is always below the pAu curve (for a given 
energy). This is because of the different flavor structure of deuterium and the fact that 
d-quarks are weighted with a factor = 1/9 in (p^, compared to the factor = 4/9 
for M-quarks. For the RHIC energy of y/s = 200 GeV, we calculate only for xp > 0.5 to 
make sure that the fluctuation lifetime significantly exceeds the nuclear radius. At the very 
high LHC energy of ^/s = 5.5 TeV, the coherence time is much larger than the nuclear 
radius for any value of xp (except at the very endpoints). Thus, the entire Xi^-range is 
shadowed. Shadowing is especially strong at LHC energies at large xp, where X2 can become 
as small as X2 ~ 10~^. At such low X2, the effects of gluon shadowing leads to a sizable 
additional suppression of the DY cross section. Without the gluon shadowing contribution 
in (p!3D, shadowing of the DY reaction at LHC would be strongly underestimated. The mass 
dependence of shadowing in DY is shown in the plot on the right-hand side of fig. ^ The 
weak dependence on M reflects the fact that shadowing for DY is a leading-twist effect, just 
as for DIS. Indeed, configurations with a — s> 1 in (|l]) and (|l^) are the analog of Bjorken's 
aligned jet configurations in DIS [0, which make shadowing persist as M ^ 00. 

We also investigate the A- and the impact parameter dependence of shadowing, with the 
results shown in fig. The amount of shadowing, i.e. the difference from unity in fig. ^ 
is to a good approximation proportional to A^^^. The deuterium curves in fig. |^ do not of 
course go to unity at A = 1 or 6 ^ 00; the flavor suppression remains in these limiting cases. 
We point out that it is a special advantage of the dipole approach that it naturally predicts 
the impact parameter dependence of nuclear effects. To obtain the b dependence, one simply 
eliminates the b integral in (|TB|) and divides by the nuclear thickness T{b) instead of A. 



5 Nuclear modification of the DY transverse momen- 
tum distribution in pA- and X)A-collisions 

The differential DY cross section as a function of the dilepton transverse momentum qp can 
be calculated in the dipole formulation as well. At the energies relevant for RHIC and LHC, 
the transverse momentum distribution of DY pairs from p{D) + A collisions can be written 
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Figure 8: Shadowing for the total DY cross section as a function of A^^^ 
and impact parameter b. Both figures are calculated for M = 4.5 GeV 
and xp = 0.5. In each plot, the lower pair of curves is for deuterium 
- gold scattering (DAu) and the two upper curves are for proton - gold 
collisions (pAu). 



in frozen approximation [|T^, p!3|] , 



Xi 



d^a{pA IIX) 
dM^ dxp d^qr ~ SttM^ Xi + 



\ da{qA 7*X) 



dlnad'^qT 



(16) 

in analogy to Eq. (|T5|) . The differential cross section for heavy-photon radiation in a quark- 
nucleus collision was derived in [O, 



da{qA 7*X) 
d In ad^qT 



1 



(27r)2 
1 



d'^pid'^P2 exp[igT ■ (pi - P2)]%*„{a, pi)'^^,g{a, P2) 



a Japi,X2) + a Jap2,X2) - o" -(a|pi - p2\,X2) 



(17) 



In this expression, the 2;-axis is parallel to the momentum of the projectile quark. Integrating 
Eq. ([1^) over q^ we arrive at the cross section Eq. (|lD. Three of the four integrations in ([T7| ) 
can be performed analytically for an arbitrary form of the dipole cross section a^- The 
details of calculations are moved to [Appendix Ij| . Since the remaining integration still has 
to be performed over an oscillating function, the gr-range in which numerical calculations 
are feasible is limited. We calculate up to qr = 10 GeV, which covers the experimentally 
interesting region. 

As in the preceding section, we perform calculations for pA and for DA scattering. Our 
results for the differential DY cross section for transverse and longitudinal pairs are shown 
in fig. 1^. We show curves only for pA in this figure, because the difference between pA and 
DA is hardly visible on the logarithmic scale. Even though X2 increases as function of qx, 
at qx = 10 GeV the coherence time tc is still larger than the nuclear radius. As already 
observed in ||13| , the differential cross section does not diverge at zero transverse momentum 
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Figure 9: Absolute value of the differential DY cross section for pAu 
scattering, divided by A{= 197). The curves are predictions for RHIC 
(^/s = 200 GeV) and LHC = 5.5 TeV) for the dilepton mass of 
M = 4.5 GeV and xp = 0.5. The upper two curves show the sum 
of transverse and longitudinal cross section (T + L) and the lower two 
curves (L ) separately show the longitudinal cross section. 



because of the flattening of the dipole cross section at large separations. On the partonic 
level, we reproduce the same asymptotic behavior that is expected in the standard parton 
model, 

darigp 1*X) 1 

oc^ forg^^oo (18) 

and 

daiiqp ^ 1*X) 1 

rfV ""4 ^^'^ 

Embedding the partonic cross section (|17D into the hadronic environment as in (|16]) will lead 
to a somewhat steeper decay at large qx, reflecting the decrease of the structure function 
F2 of the projectile with xi and the increase of Xi with qt- While this effect is quite strong 
at xp = and fixed target energies, it becomes negligible at large xp and the high cm. 
energies at RHIC and LHC. Indeed, even at qt = 10 GeV, the asymptotic limit is not yet 



fully reached. The qt dependence is still slightly less steep than in (|T8D 



To see the effect of nuclear shadowing and antishadowing we divide the nuclear differential 
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cross section by A times that of the nucleoli (2 A for DA scattering). Then, nuclear effects 
manifest themselves as a deviation from unity. The results of calculations (see [Appendix Bp 



for gold at the energies of RHIC and LHC are presented in fig. |T0| for the unpolarized (top) 
and longitudinally polarized (bottom) DY cross section ratios. Also, the difference between 
pA (left) and DA (right) now becomes clearly visible. As already explained in Sect. |, this 
difference is due to the larger abundance of (i-quarks in deuterium. Note that we neglected 
nuclear effects in deuterium and assumed isospin symmetry for the DA curves. 
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Figure 10: Nuclear ef- 
fects on the DY trans- 
verse momentum distribu- 
tion. Curves show the 
DY cross sections for pAu 
(left) and deuterium - gold 
(right) collisions per nucleon 
divided by the DY cross 
section from pp scattering. 
Solid curves are predictions 
for RHIC (.fs = 200 GeV) 
and dashed for LHC ( y/s = 
5.5 TeV) Calculations are 
for the same kinematics as 

fig- a 



For low transverse momenta, we expect DY dilepton production to be shadowed. Note 
that shadowing for longitudinal qq pairs is smaller than for transverse pairs because the lon- 
gitudinal cross section is dominated by small distances in the dipole cross section. However, 
gluon shadowing, whose onset we observe in the RHIC predictions and which should become 
the dominant effect at LHC, is about the same for longitudinal and transverse pairs. Indeed, 
we predict rather different shadowing effects for longitudinal and transverse dileptons at 
RHIC, but about the same at LHC. 

It is interesting that the effect of antishadowing, the so-called Cronin enhancement pre- 
dicted in [|l^ , disappears at the energy of RHIC after inclusion of gluon shadowing, which 



was disregarded in [T^. This reminds one of the missing Cronin enhancement in charged 



particle multiplicities that was measured at RHIC |43|. However, the RHIC data cannot be 



explained by gluon shadowing, because the x of these data is too large. Some antishadowing 
is still possible at large ~ 10 GeV at the energy of LHC as a result of the substantial 
rise of the dipole cross section with energy and the corresponding relative enhancement of 
the multiple interactions responsible for the Cronin effect. This expectation is confirmed 
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by fig. |rT|, which shows the results of calculations with (solid curves) and without (dashed 
curves) gluon shadowing. 
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Figure 11: The influence of 
gluon shadowing on the DY 
cross section. Dashed curves 
are calculated without gluon 
shadowing, i.e. Rq = 1 in 
(\r^), while solid curves in- 
clude gluon shadowing. The 
influence on the longitudinal 
DY cross section is shown 
separately in the two left- 
hand plots {L). The two 
plots on the right show the 
DY ratio for the sum of the 
transverse and longitudinal 
cross sections {T + L). Cal- 
culations are for the same 
kinematics as fig. ||. 
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The DY process with the production of a longitudinally polarized photon manifests 
stronger effects of antishadowing (fig. |10| bottom and fig. |Tl| right), as was earlier observed 
in 1 12]. However, this enhancement of the longitudinal cross section will be difficult to see 
in experiments because the transverse cross section is so much larger than the longitudinal 
one. All nuclear effects are expected to vanish at very large qt- 

One can also study the moments of the transverse momentum distribution. A frequently 
measured characteristic of nuclear effects is the broadening of the mean value of the DY 
transverse momentum squared, which is the difference between the values of mean transverse 
momentum squared measured in pA and pp collisions, 

HQt) = (4) A - {ql)N , (20) 

where 

da^/d^qx is the proton-nucleus DY cross section given by Eq. |16[ 

It is easy to understand from fig. |IO|that a nuclear target leads to a larger mean transverse 
momentum of DY dileptons than a proton target: low qt pairs, corresponding to large 
arguments of the dipole cross section (ap), are shadowed, while high transverse momentum 
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pairs are not (color filtering [2C]). However, the actual numerical value of broadening, i.e. the 
increase of the square mean transverse momentum, depends on the maximum included in 
the analysis. This is not an artifact of our approach, as this is also the case in experiment. 

According to (|1^) the numerator in (|21| ) diverges at large qx for the transverse cross 
section. Even after averaging over the projectile parton distribution, the integral is very 
slowly converging, and one has to introduce an upper cutoff g^"^ since there is a maximal 
transverse momentum accessible in experiment. On the other hand, the large-g^ tail of the 
differential cross section should be the same for nuclear and nucleon targets since no nuclear 
effects are expected at large gy- For this reason one may think that the divergence would 
cancel in the difference in Eq. (|20|) and render the result cutoff independent. This might be 
true if no nuclear effects occurred in the integrated DY cross section, i.e. if a^j^ = Aa^. 
However this is never the case. For long coherence time, tc ^ Ra, shadowing diminishes 
the DY nuclear cross section, i.e. the denominator in the first term in Eq. pPj). As a result, 
the high-gr tail of the nuclear distribution is renormalized and undercompensated by the 
second term in Eq. (]20|). This is why there is sensitivity to the upper cutoff q^"'^ in our 
results, and this sensitivity is even more pronounced at the higher energies where shadowing 
increases. On the other hand, for short tc ^ 0, where shadowing vanishes, energy loss has a 
similar effect of suppressing the DY cross section on nuclei [^, |^, |^ . 



Note that in at least some experiments P, ^ the transverse momentum broadening 

(22) 



is extracted from the data by fitting the points with the functional form 



(Pqj 



where typically n = 6 for both a proton and a nuclear target. The fast decay of the differential 
cross section in Eq. (^21) is due to the increase of Xi as function of g^- As mentioned above, 
this effect vanishes at large xp and high s. The gy-distributions at RHIC and LHC (at large 
xp) will be much flatter than (|2^). 

In fig. O we compare A-dependences of the broadening 6{q^) calculated with different 



cutoffs, qp'^^ = 5 (bottom curves) and 10 GeV (upper curves). The main observations are 
first that broadening is roughly proportional to the path length in the nuclear medium. This 
is true for any value of g^°^. Furthermore, 5{q^) can become quite large for heavy nuclei, 
around 1 GeV at RHIC and around 3 GeV at LHC. While the influence of gluon shadowing 
on broadening is rather weak, the g^"^-dependence is quite strong. This is studied in more 



detail in fig. |13[ Increasing the transverse momentum cutoff from 5 GeV to 10 GeV at 
RHIC energy leads to an increase of nuclear broadening of slightly more than 50%. At LHC 
energies, however, where one still has nuclear effects in the transverse momentum distribution 
at rather large values of qp, broadening increases by a factor of 3. Therefore, the DY process 
turns out to be a less than ideal tool to measure the broadening of the transverse momentum 
distribution for a quark propagating through nuclear matter. 

Finally, we calculate the energy dependence of nuclear broadening, shown in the plot 



on the right in fig. [T^. Again, calculations are performed for two different values of the 



transverse momentum cutoff. Note that the shape of the curve depends strongly on g^' 
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Figure 12: Nuclear broadening for DY dileptons from pAu at RHIC and 
LHC. Broadening depends on the transverse momentum cutoff qjp"'^ . In 
each plot, the lower pair of curves is calculated for g^"^ = 5 GeV, while 
the upper pair is for q^"-^ = 10 GeV. The solid curves include gluon shad- 
owing, the dashed ones do not. Calculations are for the same kinematics 
as fig. |. 



For a transverse momentum cutoff of 5 GeV, there is almost no energy dependence of S{q^) 
above RHIC energy. The situation looks different if the transverse momentum cutoff is 
10 GeV. In this case, broadening does increase as function of energy. It will therefore be 
difficult to draw conclusions from possible future data on the energy dependence of S{q^), 
since presumably one will only see the cutoff dependence. 

Valuable insight into the relationship between shadowing and broadening is gained if one 
performs the integration in the numerator of ( ^T]) analytically for q^°'^ oo. Without the 
projectile parton distribution, the result reads 



1t)a{n) -V^ + -^a{n) 



'DY 



a 



where the wave function squared has the form given in Eqs. ( |B.2|) or 



a^m^; a'j^{a) = J d'^qx ^^^{qT, a); and 



(jA{N)^ 



X2 



(23) 

|);r/2 = (l-a)M2 + 
(24) 



Note that without gluon shadowing, all the difference between a nucleus and a nucleon occurs 
in the denominator of the second term, (T^^'*(a). Including gluon shadowing, one has 



C^(a;2) =C^^(a;2) / d'h RG{x2,h)T{h) = {x2 



A - / ARG{x2,b)T{b) 



(25) 



At the same time, the problem of the divergence at large g-j- is moved to the integral of the 
LC wave function squared in the second term, which has logarithmic singularity at p ^ 0. 
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Figure 13: The cutoff (left) and energy (right) dependence of transverse 
momentum broadening for pAu scattering. The calculations includes 
gluon shadowing, and are for the same kinematics as fig. 



Thus, the broadening Eq. (|20| ) takes the form 



I^ARdy-jJ d% ARG{x,,b)T{b)^ / , (26) 



where ARdy and ARg are the nuclear suppressions for DY and for gluons, i.e. the difference 
from unity, ARdy = 1 — cy'j^ / {Aa]^) (see Eq. |]for ARg). 

We reach the interesting conclusion that if there is no shadowing, the broadening of the 
transverse momentum vanishes as well. This is a manifestation of a close relationship between 
broadening and shadowing in the regime of tc ^ Ra- Indeed, broadening is interpreted 
in the LC dipole approach as color filtering. Namely, the mean size of a qq dipole that 
can be propagated through nuclear matter decreases with the size of the configuration; 
correspondingly, the intrinsic transverse momentum of the surviving configurations rises. 
Shadowing occurs due to the same phenomenon. Once one says that shadowing is negligibly 
small, this also means that the dipole is too small to undergo multiple interactions. However, 
in this case no color filtering occurs either. This rather intuitive result looks very nontrivial 
in the framework of parton model, where one can get broadening even without shadowing 



Note that gluon shadowing seems to reduce the amount of broadening in Eq. ( |26| ) . This 
should be expected, since gluon shadowing reduces the nuclear thickness, cf. Eq. ([T3|) , and a 
more dilute medium leads to less broadening. However, shadowing for DY, ARdy, increases 
with gluon shadowing. Numerically, we find that the influence of gluon shadowing cancels 
to a large extent and that broadening is almost independent of Rg, see Fig. |T^. 

In Eq. ( pTD we avoided the (logarithmic) divergence in (q^), related to the singular 
behavior of Ki(?7p) at small rjp for transverse photons, by introducing above an upper cutoff 
q^p on the integrals over qx- These numerical results constitute quantitative predictions of 
the LC target rest frame formulation that may be compared to experiment noting that g™"^ 
is a physical parameter related to the acceptance of the spectrometer in the measurement. 
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Motivated by the desire to understand these same numerical results analytically, we next 
examine the theory employing certain simplifications and approximations. 

Our main approximation is to make the replacement of the fiuctuation distribution by a 
Gaussian, 

- + (1 - afWn\-^l^\ (27) 

where = 2(3iff, with and (3 chosen to give an acceptable match to the actual fiuctuation 
distribution, which we simplify to be 

P)? = + (1 - afWKl{^p), (28) 

noting that is small and the longitudinal contribution is about a 10% correction to the 
momentum distribution. Because the distribution of fiuctuations is integrable with the cutoff, 
the average size of the fiuctuation is meaningful, and with the Gaussian wave function, the 
mean-square transverse spread of the fiuctuation is (p^) = 1/A;q. 

With such approximations, the integral in Eq. ( p!7D can be evaluated analytically. The 
result of carrying out the integrals over pi, p2, and qr gives 

2 o d^afy ^ /• ,2, \2C{a)a'TA{h) , 2 [C{a)a^TA 
" Q.T Qt "771 7779 — ~ TT / a 



d{\na)(PqT 7 [ [2 fc^ + C(«) «2 7^a(&)]' 

,/,.6^!WiiM + M^™!)]!+,,,y. (29) 

/cq 2 Kq 



where we have omitted the prefactors in Eq. (|27|) , and the expansion is useful for examining 
the limit of weak shadowing. Likewise, 

'^d{lna)d'qT kl J 2kl + C(a) a^T^ib) 

J y 2 /cq J 

In these expressions we have introduced an effective C = C{a), defined independent of (3 so 
that the Drell-Yan cross section for the GW color dipole cross section is reproduced in the 
approximation to it on a nucleon, 

C{a) = I^}M/}-P)d'P, (31) 
J e oP a^p^d^p 

We could optimally match approximate theory to the exact one by choosing and (3 
numerically for a given g™"-^, but since our interest is insight rather than numerical precision 
at this point our conditions are simply the following: (1) Preserve the integral j p'^'K\{rip)d^ p, 



OO /"OO 

2„2 



p'Kl{7^p)d^p = / p'e-'^oP'd'p, (32) 
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determining = 16/3^/3. (2) Adjust /3 to preserve (p^) using the asymptotic form for 
Ki(a;). This gives (3=1. We determine the effective momentum cutoff corresponding to 
this Gaussian by comparing the exact numerical value for tv for the nucleon to the 
approximate value. The latter is obtained from the small- A limit of Eqs. (pQ]) and (pO|), 



2l)N = 2{a\l + (1 - af))/{a\l + (1 - af)/kl), 



(33) 



where the brackets indicate the convolution (16) with the quark distribution function. We 
find that for g^"^ = 10 GeV, {qT)N agrees with the RHIC values to about 10% and LHC 
values to about 30%. 

For a nucleus, the integrals over impact parameter in Eq. (|29| ) and Eq. (|30| ) may be 
carried out analytically for a sharp-surface density model in which the density is constant 
at po out to radius R1/2, and zero beyond. Then, from Eq. (^), we find 



/ 



j2 2 

d qtQt 



DY 



d{\n a) d^qT 



(34) 



and from Eq. (^ 



/ 



d^QT 



d^a 



qN 
DY 



/2 



d(\iaa)d'^qT l/^^o 



ln{l + y) + — - y 



(35) 



where y = C {a)a'^ poRi/2/ k^. To obtain (q^) for the nucleus, it is necessary to perform a 
convolution of these expressions with the quark distribution function of the nucleon projectile 
as in Eq. (|I6|). 

It is interesting to examine the expansion of Eq. (0) and Eq. ( ^5]) in powers of y. 
Recalling that /cq = l/(p^), and further noting that L = 2i?i/2 is the diameter of the 
nucleus, this is essentially an expansion in y = L/2X where A = l/{aqq)po is the mean-free 
path of the fluctuation. For most values of a, the number of interaction mean-free paths in 
crossing the nucleus is tiny due to the large value of in 77^ = (1 — a)M^ -|- a^m^. However, 
for a ^ 1 , the mean transverse separation of the fluctuation may become relatively large (we 
find (p^)^/^ = 0.7 fm for rriq = 0.2 GeV). In fact, for a large nucleus (y4=208), we see that 
for RHIC energies C(q! ~ 1) = 2.5 and y ^ 1.5. Clearly, these larger fluctuations have a 
relatively small mean-free path A and are subject to appreciable color filtering in traversing 
the nucleus. For LHC, C(a ~ 1) = 5.7 and y ~ 3.4. At he same time, the amount by which 
(q^) differs its value for a nucleon, Eq. (|33|), grows and hence nuclear broadening also grows. 
In this fashion, the physics of nuclear broadening is again seen to be directly related to color 
filtering for a ^ 1. Since the expansion of Eqs. ( P^ and ( |55| ) converges slowly in the region 
where the largest contributions to the nuclear broadening occur, it is necessary to evaluate 
the integrals over b and a without making an expansion to calculate S{q^) with sufficient 
accuracy. 

With these approximations (we have also omitted gluon shadowing, which has been 
shown to have a weak effect on 6{q^)), we find for heavy nuclei tc ^ Ra, S{qT) ^ (^1); 
i.e. S{q^) is very nearly linear with A^^^, just as we found in our more exact numerical 
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studies. The constants of proportionality are about 2.2C(x2) at RHIC and l.lC{x2) at 
LHC, which overestimate the slope of the exact results by about 20% and 44%, respectively. 
Although linearity in A^^^ would follows in the weak scattering limit, we again remark that 
our analytical calculations indicate substantial effects from higher-order multiple scattering. 

In the opposite limit of short tc ^ 0, broadening of a quark is known to rise linearly with 
the length of the path of the quark in nuclear matter before the DY reaction occurs [^9|, |l^ , 



(36) 



where (Ta) = J (PhT\{b) / A is the mean nuclear thickness. At the same time, broadening 
for the DY pair, which carries only fraction a of the quark momentum, is reduced. 



HQt)dy 



where 



1 ^ 
/q da o? J (Pp ^q^* (p, a, M^) cr^-(p, x) 



Jo da J (Pp 



^„..(p,a,M2) <-(p,x) 
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(37) 



(38) 



This value depends on M logarithmically, and for M ~ 5 GeV, (a^) ~ 0.9. 

It is interesting to compare the nuclear dependences of broadening in the two limiting 
regimes tc 3> Ra, Eq- (PD, and tc — ^ 0, Eq. (^61). Expanding the nuclear cross section 



Eq. (|^) in a^^TA^h), we can then perform the integration for longitudinal DY photons (for 
transverse photons it diverges). Then we arrive at the same expression Eq. (^), except tat 
it acquires an extra factor 



K 



2 ' 



where 



(39) 



(40) 



Applying the p^ approximation Eq. (|ToD for the dipole cross we get K = 4/5. Thus, the 
broadening for longitudinal DY pairs in the asymptotic regime tc ^ Ra matches rather well 
the low energy regime. 

The dependence of the broadening S{q'j^) on the cut-off is an unpleasant property that 
brings uncertainty to the comparison of theory with data. As we mentioned, it is related to 
the large qt behavior Eq. ( [T^ ) of the DY cross section, leading to a logarithmic divergence 
in the integral over weighted with q'^, 



pA 
DY 



<1t)= j d^qTqla^oyilT^a) , 



(41) 



which is the numerator in ([2T|). Since this integral has exactly the same divergence at large 
^max nuclear and nucleon targets, it must cancel in the difference. 



5{aDY qr) = W^dy It) - ^ Wdy It) 



(42) 
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and the result should be independent of the upper cut-off qjp"-^ when it is sufficiently large. 
One can also normalize this difference dividing both terms by A {(Tdy)- Unfortunately, the 
result is not an exact measure of the broadening of the transverse momentum of a quark 
propagating through a nucleus. However, these quantities are independent of the experi- 
mental acceptance and this fact makes it a better observable than the broadening, 
Eq. (^0|), to compare with theory. 



6 Polarization of DY pairs 

In the preceding section, we separately calculated the DY cross section for transverse and 
longitudinal photons. In experiment, the different polarizations can be distinguished by 
investigating the angular distribution of DY pairs. The most general form of the DY angular 



distribution reads 21 



, ,,.2 J n,. « 1 + ^ cos2 e + ^ sin(2e) cos(0) + sm\e) cos(2</.), (43) 
ax F d M ^ a cos 6 acp 2 

where 9 is the angle between the muon and the 2;-axis in the rest frame of the virtual photon, 
and (j) is the azimuthal angle. Of course, A and (f) depend on the choice of 2;-axis in the dilepton 
center-of-mass frame. Since the dipole approach is formulated in the target rest frame, it 
is convenient to put the 2;-axis in the direction of the radiated photon |TT|. The target rest 



frame and the dilepton center-of-mass frame are then related by a boost in ^-direction, so 
that the transverse polarizations are the same in the target rest frame and in the photon 
rest frame. Note that in the dilepton center-of-mass frame, the 2;-axis is antiparallel to the 
target momentum. This frame is called the u-channel frame and the curves we present are 
valid in this frame. 

The 0-dependence of the cross section is difficult to measure. At RHIC, only the value 
of A can be measured Since A = for transverse and A = —1 for longitudinal, one 



obtains after integration over the azimuthal angle. 



A = (44) 



Our results for A as function of qx at RHIC and LHC energies are shown in fig. Note 
that = X1X2S — is a Lorentz invariant variable and therefore has the same value for 
all choices of the 2;-axis. At low g^, A shows a clear deviation from unity. This deviation 
increases with energy, because the longitudinal cross section is more strongly dominated by 
small distances than the transverse. Since the dipole cross section grows faster with energy 
at small separations, the relative amount of longitudinal photons at small qx increases with 
energy, see fig. ^. At very large transverse momentum, A will eventually approach unity 
because the longitudinal cross section falls off steeper than the transverse one (p!8|). It is 
also interesting to investigate nuclear effects on A. For pp collisions, A has been calculated 
previously in [T^. One sees from fig. that the nucleus leads to a stronger suppression 



at small q^ than a proton target. This is because longitudinal photons are less shadowed 
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Figure 14: The parameter A in DAu collisions at RHIC and LHC 
(left). The figure on the right illustrates nuclear effects on X, by showing 
the ratio of A for DAu scattering and for pp scattering. Calculations are 
for the same kinematics as in fig. |^ and include gluon shadowing. 



at g-r — ^ than transverse, fig. |T^. Nuclear effects on A vanish at very large transverse 
momentum, as expected. 

An interesting feature of our result is that the parton model relation 



l-A-2i^ = 0, 



(45) 



which is known as Lam- Tung relation [|2T|, is not fulfilled in the dipole approach. The 
violation of the Lam- Tung relation becomes apparent from the behavior of A as ~^ 0. 
Kinematics require that the double-spin-flip amplitude u vanishes at qt = 0, thus A should 
approach unity. Remarkably, this relation holds to order as in the QCD improved parton 
model. Note that ( ^Sj) even holds separately for the annihilation and for the Compton 
corrections in the parton model. Even at order a^, the corrections to (|45|) are very small 
PI , by order of magnitude 1 — A — 2z/ ~ 10~^. Experimentally, the Lam- Tung relation is 
clearly violated ||5T|, 1 — A — 2z/~ —1. However, all experimental data lie at large X2 and 
are therefore not relevant for our calculation. At present, the reason for the violation of 
(^) is not known. Possibly, higher twist effects, included by default in our dipole approach, 
can provide an explanation [Q. However, it is not yet possible to apply the resummation 
technique |Q for logs in (g-r/M) to higher twists, which would be necessary to compare to 
future data on DY angular distributions from PHENIX. 

The dipole formulation provides a much easier way to calculate the DY transverse mo- 
mentum distribution even at low qt- It is not surprising that the result from the 0{as) 
parton model calculation is not reproduced in our approach, because the dipole picture is 
not an expansion in powers of the strong coupling constant. Instead, all contributions from 
higher order graphs that are enhanced by a large log(l/x2) are contained in the dipole cross 
section. Moreover, if we use a phenomenological parameterization of cr^, some higher twists 
and nonperturbative effects are contained as well. The Lam- Tung relation is violated in our 
calculation because of the flattening of the dipole cross section. Indeed, for a quadratically 
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rising cr^g{p,x) = C{x)p^, A would vanish at qr = 0, as can be seen from (0). Note that 
the remaining coefficients for the DY angular distribution in the dipole formulation 



can be obtained from Eqs. (22)- (25) of |TT|]. It is easy to check that (^) is fulffiled only 
with cr^(p, x) = C{x)p^. Thus, the Lam- Tung relation is also violated in pp-colhsions [ p!3[ . 
Finally, we would like to stress that the behavior of A at low qx depends heavily on the large 
p behavior of the dipole cross section, which is not well constrained by DIS and diffractive 
data. However, even if (|^) is unrealistic at large separations, the qualitative behavior of A 
will remain the same for any flattening parameterization. We therefore believe that it is 
worthwhile to measure A at RHIC, because such data could give us information about the 
dynamics beyond the conventional parton model. 



i^^f(^i,^2)^ .?y;;L ^^ > (46) 



7 DY process in heavy ion collisions 

As we pointed out in the introduction it is impossible to predict nuclear shadowing for the DY 
process in nuclear collisions from the parton model using data for DIS and the DY reaction 
on proton and nuclear targets. Indeed, the nuclear effects to be predicted are presented in 
the form, 

ABa^^{x^,X2) 
where 

a^f(xi,X2) = nY,Z} I d% I dhTA{s)TB{b - s)[q}{x,)R^{x,,s) q}ix2)Rf{x2,b - s) 
J J J 

+ q}{xi)Rf{xi, s) q}{x2)R^{x2,b- s) + q}{xi)Rf{xi, s) q}{x2)Rf{x2,b - s) 
+ q}{x,)Rf{x^,s)q}{x2)Rf{x2,b-s)'\ , (47) 

a^^{xuX2) = Nj2zj[q}{xi)q}{x2) + q}{xi)q}{x2) 

f 

+ q}{xi) q}{x2) + q}{xi) q}{x2)\ , (48) 

and b is the impact parameter of the beam [A) and target [B) nuclei. The ffist two terms in 
Eq. (^Vp correspond to (i) annihilation of a valence quark of the beam with a sea antiquark 
of the target; (ii) a sea antiquark of the beam with a valence quark of the target. The third 
and fourth terms correspond to (iii) both the quark and antiquark originating from either 
the sea of the beam or target. Provided that the Drell-Yan /^-factor is independent of A 
and B, the overall normalization factor in these expressions is irrelevant for shadowing 
Eq. ( 1^ ) since it cancels in the ratio. This assumption is supported by experimental data on 
the K factor for various heavy ion collisions ||54|| . 

Note that a generalization of Eqs. (|I6|) - (p^ to the case of AB scattering including higher 
twist effects, would be much more complicated than Eqs. (0) - (|48|). Such higher twist 



effects, however, will be important only for nuclear effects in the DY transverse momentum 
distribution in AB, which are beyond the scope of this paper. 
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All parton distributions in Eqs. (^) - (^8]) are taken at the same virtuality = M^. 
For the sake of simplicity, we assume that the shadowing, Rf, is the same for sea quarks and 
antiquarks, and we neglect the isospin noninvariance of the sea distribution at moderately 
small X ~ 0.1 p^, which can easily be taken into account. 



The nuclear shadowing Rq-^{x, s) is given in Sect. |7.1| below. Note that -R^(x, b) at small 



X is a function of nuclear thickness. It vanishes at large impact parameters on the nuclear 
periphery, but it reaches its maximum at 6 = 0. Data for DIS or the DY reaction in 
pA collisions provide information only about the 6- integrated shadowing effect. Knowledge 
of only such integrated shadowing is not sufficient for the calculation of shadowing in an 
AB collision, Eq. (|47D . Nevertheless, it was assumed in ^ that shadowing is independent 
of impact parameter, R^{x,Q'^,b) = R^{x,Q'^). Clearly, such an ad hoc assumption is not 
justified and leads to basic consequences which cannot be accepted. For instance, it precludes 
any dependence of shadowing effects on centrality in heavy ion collision. 

As we demonstrated above, the LC dipole approach provides direct access to the im- 
pact parameter dependence of shadowing effects. In the following subsection, we calculate 
shadowing for valence and sea quarks and compare the 6-integrated result to the EKS98 
parameterization Of course, the impact parameter dependence of shadowing is taken 



into account in Sect. 7^, where we predict shadowing for the DY reaction in a nuclear 
collision in the integrated form Eq. (0), as well as a function of centrality. In order to 
calculate -R^y(xi, X2, &) for a collision with impact parameter b, one should just eliminate 
the integration over b in the numerator of Eq. (^), and replace 



AB^TABib) = / d^sTAis)TBib 



(49) 



in the denominator of Eq. (|^). Comparison with the minimal-bias events for DY dilepton 
production in heavy ion collisions would serve as a rigorous test of the theory. 



7.1 Nuclear shadowing for sea and valence quarks 



Since we rely on the factorization relations Eqs. (|46|) - (^8]), we can calculate shadowing 
Rf{x,Q^) and R^{x,Q'^) in DIS, which looks somewhat simpler as it does not include a 
convolution with the initial quark distribution. 

Shadowing for sea quarks is calculated with Eq. ([T3|) and is given by 



Rs{x,Q\b) 



(50) 

The light-cone wavefunctions \E'qg for the transition 7* qq can be found in the literature, 
see e.g. [|T7|. Let us recall that this expression is valid only for the so-called "frozen" ap- 
proximation, i.e. in the asymptotic regime tc ^ Ra, which takes place at very small x. In 
the transition region < Ra one should employ the LC Green function technique. This 
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(1-2^ ^)Rg{x, ^Ip'^ b)TA{b)f 


TA{b) da J d^p 
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was done in [5(:, 39|, although gluon shadowing was neglected in those calculations. At very 
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small X gluon shadowing is essential, as demonstrated above (see comparison with data from 
the NMC experiment in [0). 

Nuclear shadowing for valence quarks has never been calculated. This shadowing is 
usually believed to be small 0], if it occurs at all. We shall demonstrate, however, that 
shadowing for valence quarks is quite sizable, even stronger than the shadowing of sea quarks. 
This is another new result of the present paper. We note in this regard that the nuclear 
structure function F2{x) is different from the quark distribution function in an essential 
way; namely, the former contains shadowing effects and therefore the baryon number sum 
rule is not applicable to it p4|, a difference that may explain the discrepancy of our results 



compared to Ref. 

Note that we call the ratio Eq. (pOf) shadowing for sea quarks because the dipole cross 



section crg^(p, x) includes only the part that rises with energy, corresponding to gluonic 
exchanges in the cross-channel. Therefore, this is the part of the sea generated via gluons 
(there are also other sources of the sea, for instance the meson cloud of the nucleon, but 
they vanish linearly with x or faster). The fact that the color-dipole cross section includes 
only the part generated by gluons is the reason why it should be used only at very small 
X < 0.01, where the sea dominates. This part of the dipole cross section can be called the 
Pomeron in the language of Regge phenomenology. In the same framework, one can relate 
the valence quark distribution in the proton to the Reggeon part of the dipole cross section, 
which has been neglected so far. 

Thus, to include valence quarks in the dipole formulation of DIS, one should replace 

a^^ip, x) => af^ip, x) + a^ip, x) , (51) 

where the first (Pomeron) term corresponds to the gluonic part of the cross section, which 
we have used so far. It is responsible for the sea quark part of the nucleon structure function. 
The second (Reggeon) term must reproduce the distribution of valence quarks in the nucleon; 
this condition constraints its behavior at small x. One can guess that it has the following 
form, 

a^ip,x)=Np^^ , (52) 

where ^/x should reproduce the known x dependence of valence quark distribution (as, in 
fact, motivated by Regge phenomenology), and the factor is needed to respect the Bjorken 
scaling. The factor N will cancel in what follows. 

We are now in a position to calculate shadowing for valence quarks by inserting the cross 
section Eq. ( pT] ) into the eikonal expression Eq. (0). Clearly, only the gluonic part of the 
dipole cross section is subject to gluon shadowing, i.e. only the first (Pomeron) term of 
Eq. (^) should be multiplied by Rq- Furthermore, if one expands the numerator in powers 
of and picks out the linear termQ, then one arrives at the following expression for nuclear 

^The small size of a^{p, x) at small x motivates such an expansion; however, one should note that it would 
not be proper to include the higher powers of the Reggeon cross section. Indeed, the Reggcons correspond to 
planar graphs. These cannot be eikonalized since they lead to the so-called AFS (Amati-Fubini-Stangelini) 
planar graphs, which vanish at high energies pS]. 
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Figure 15: Shadowing for sea 
and valence u-quarks in DIS off 
gold at Q = 4.5 GeV. Solid lines 
are calculated from Eqs. ([5^ - 
(|5^, while dashed curves show 
the EKS98 parameterization 



shadowing of the valence quarks, 
R,{x,Q\h) 
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(53) 

This shadowing is even stronger than the shadowing for sea quarks, Eq. (^Up. Indeed, for 
weak shadowing we can also expand Eqs. ( ^0]) and (|53|) in powers of a^. Then one obtains 
a shadowing correction 1 — = ^aeffT^ in (^) which is twice as large as that for sea 
quarks I - Rs = ja^jfTA in (0). Here a^fj = {a'^g) / {crqq) 0. 

These estimates rely, however, on the p^-approximation for the dipole cross section and 
on the assumption that shadowing is weak. The result of our calculation, including gluon 
shadowing and a realistic parameterization of the dipole cross section, is shown in fig. ^ 
We show only the 6-integrated shadowing, which is given by 
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for valence quarks and by 
Rs{x,Q^) = 



(55) 
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for sea quarks. Shadowing for valence quarks is still stronger than it is for sea quarks, but not 
by a factor of 2. However, valence quark shadowing calculated in the LC approach is much 
stronger than in the parameterization of [Q. Note that the authors of ^ force their fit to the 
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nuclear valence quark distribution to fulfill the baryon number sum rule, which is however 
violated by shadowing effects. Therefore valence quark shadowing is underestimated in 

Unfortunately, it will be impossible to extract the \ow-x valence quark distribution of 
a nucleus from DY experiments, because the nuclear structure function is dominated by 
sea quarks. Maybe neutrino-nucleus scattering experiments could provide data that help to 
differentiate shadowing for sea and valence quarks. 



7.2 Modification of the DY cross section in heavy-ion collisions 



We can now make use of Eqs. ( [46D - (|48| ) and predict nuclear effects for the cross section of 
DY lepton pair production in heavy ion collisions. We perform calculations at large xf, where 
the structure function of the target nucleus, say nucleus B in (0), enters at small 0:2 -C 0.1 
and is therefore subject to shadowing. We calculate the shadowing ratios Rs^v in i W^ 
function of impact parameter from Eqs. (|50D and (0). However, shadowing considered in 
the previous section is not the only nuclear effect affecting the ratios X2). Indeed, 

the parton distributions of the projectile nucleus. A, are sampled at xi > 0.1, i.e. in the 
region where the EMC effect must be taken into account. This leads to a 10 — 20% 



suppression at medium-large Xi and a strong enhancement at Xi ^ Xp ^ 1 due to Fermi 
motion. Also, a small 2 — 3% enhancement is known to exist at Xi ~ 0.1. All these are 
medium effects caused by the difference between the properties of bound and free nucleons. 
Since the nuclear density, apart from its surface, is approximately the same for all nuclei 
except for the lightest ones, one can assume that these medium effects are about the same 
for all bound nucleons in all nuclei. This assumption is supported by data, which displays no 
strong A-dependence of the EMC effect from medium through heavy nuclei [0]. Of course, 
there is no contradiction here with our previous statement about a strong impact parameter 
dependence of nuclear shadowing. 

For our actual calculations, we employ the EKS98 0] parameterization of the nuclear 
parton distributions for nuclear effects at large x (EMC, Fermi motion). For the parton 
densities of a proton, which are needed as a baseline for EKS98 and for calculation of the 



denominator (|48|) , we use the CTEQ5L parameterization Note that, unlike in the pA 
case, parton densities of the proton are not only needed at large x, but also at very small 
X -C 0.1. The CTEQ parameterization is applicable down to x = 10~^, which is sufficient 
for RHIC {x2 ~ 0.001) but not for LHC, where values as low as x ~ 10~^ are reached. For 
the LHC calculation, we therefore switch to GRV98LO [Q, which is applicable down to 
X = 10-9. All evolution codes are taken from CERNLIB PDFLIB 8.04 gT[. 

Our predictions for the DY modification in AA collisions at RHIC and LHC are shown in 
fig. [16|. For comparison, we also display the analogous curves for pA scattering from Sect. ^ 
The plot on the left shows that DY dilepton production in AuAu collisions is suppressed 
much more strongly than in pAu collisions at the same kinematics, except for large values 
of xp, where Fermi motion makes the nuclear suppression vanish. However, the strong 
suppression in AuAu collisions is not just a combination of shadowing and the EMC effect. 
As one can see from the plot on the right, DY from AuAu is still suppressed by about 40% 
compared to pp at very large 6, where such nuclear effects are absent. Indeed, a nonnegligible 
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Figure 16: Nuclear effects on the DY process in heavy ion collisions at 
M = 4.5 GeV. For comparison, the corresponding curves for pAu scat- 
tering (fig. 0j are also displayed. The plot on the right, which shows the 
impact parameter dependence of nuclear effects, is calculated atxp = 0.5. 
The curves for AuAu collisions do not approach 1 at large b because of 
the different flavor composition of a nucleus. Note that all curves are 
divided by the values for pp scattering. 



part of the suppression in nucleus-nucleus collisions is due to the different flavor composition 
of nuclei. A heavy nucleus like gold {A = 197, Z = 79) has more neutrons than protons, thus 
the average nucleon in this nucleus has an excess of d-quarks over w-quarks, compared to a 
proton. Since (i-quarks enter the lowest-order expression for the DY cross section ( ^8]) with 
a weig ht factor of = 1/9 (compared to = 4/9) one observes an additional suppression. 



Finally, we mention that the suppression at large b could be even stronger than in fig. 16 



because neutrons are predominantly located at large impact parameter. This effect is not 
taken into account in our calculation, where all nucleons are assumed to be uncorrelated. 



8 Summary and outlook 

We have presented an analysis of nuclear effects in DY dilepton production in pA and AB 
collisions. All calculations are performed within the light cone dipole formalism for the DY 
process, because this approach suggests a very simple and intuitive treatment of nuclear 
effects. It essentially simplifies at high energies, where the coherence length for the DY 
process substantially exceeds the nuclear size. Then one can employ the eikonal formalism 
to describe multiple interactions, since the different eigenstates of interaction do not mix. 
This regime, relevant to the energies of RHIC and LHC, is considered throughout the paper. 
Formulas for the DY cross section are sufficiently simple to incorporate realistic nuclear 
densities and a realistic parameterization of the dipole cross section. The predictions for 
RHIC and LHC presented here can therefore be compared to future data and serve as a test 
of the theory. 
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Since we assume the coherence length to be long, gluon shadowing becomes important. 
The coherence length for higher Fock states is shorter than that for the lowest Fock state, 
and it is usually the order of, or shorter than, the nuclear radius. In this case, simple 
eikonalization cannot be applied to gluon shadowing and we employ the LC Green function 
formalism, which takes account of the variations in size of the projectile fluctuations as they 
propagate through the nucleus. The gluon shadowing incorporated in our calculations is 
especially important at very high energies, or small X2 <^ 1, and it makes the gg-nucleus 
cross section saturate at a value significantly smaller than the geometrical limit of 2tiR\. 

Not surprisingly, gluon shadowing leads to a stronger suppression of the total DY cross 
section than one would obtain from quark shadowing alone. Note that while fixed target 
experiments at medium-high energies find suppression of the DY cross section only at large 
Feynman-XiT', partially caused by onset of shadowing and energy loss, we expect the entire 
xp > 0-range shadowed at RHIC and LHC For the DY transverse momentum distribution, 
gluon shadowing leads to nontrivial modifications. At low lepton-pair transverse momentum, 
gluon shadowing enhances the suppression already expected from quark shadowing, but at 
intermediate transverse momentum, it strongly reduces the enhancement from the Cronin 
effect. This observation resembles the missing Cronin enhancement in charged particle mul- 



tiplicities [Q. However, the latter effect cannot be due to gluon shadowing, because the 
X of the data is too large. The observed nuclear suppression is probably caused by final 
state interactions with the produced matter and is related to the induced energy loss and 
absorption. This will be studied in more detail in a forthcoming publication. 

Furthermore, we calculate the nuclear broadening of the mean transverse momentum 
squared for DY dileptons produced in pA collisions. This quantity turns out to be divergent 
for radiation of transversely polarized DY photons. We demonstrate that this problem is 
caused by nuclear shadowing in the total DY cross section and that the two phenomena, 
broadening and shadowing, are closely related. As a result, the predicted broadening of 
the transverse momentum squared depends strongly on the upper cutoff on qx- This is not 
purely a problem of the theoretical approach, as this cutoff dependence is present also in the 
experimental analysis. We found that {5q^) can vary up to a factor of three, depending on the 
cutoff. We suggest a different observable in which the divergent tails of the qt distribution 
cancel and therefore render the result independent of the cutoff. 

We separately analyze the differential DY cross sections for transversely and longitudi- 
nally polarized pairs. In both cases, the differential cross section da/d^qT does not diverge 
as ^ 0. This result follows naturally in the dipole approach as a consequence of the 
saturating dipole cross section. In the parton model, a more complicated resummation of 
logs in q^/M is necessary in order to render the DY cross section finite as ~^ 0. On 
the partonic level and for large transverse momenta, we reproduce the behavior expected 
from perturbative QCD, namely dar/d'^qr oc q^ for transverse pairs and dai/d'^qx oc q^ for 
longitudinal. 

Experimentally, the transverse and the longitudinal DY cross sections can be distin- 
guished by investigating the angular distribution of DY pairs. We calculate the parameter 
A, Eq. (0), which characterizes the relative contribution of transverse and longitudinal pairs 
to the DY cross section as a function of the dilepton transverse momentum and, in addition. 
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investigate nuclear effects on A. Although these nuclear effects do not turn out to exceed 
6%, the different g-p-dependence of the transverse and longitudinal cross sections leads to a 
nonmonotonic behavior of A that can be observed in future experiments. 

Finally, we present estimates for nuclear effects in heavy ion collisions at the energies 
of RHIC and LHC We make use of QCD factorization, but calculate nuclear shadowing 
for sea and valence quarks separately within the LC dipole approach. Contrary to usual 
expectations, we found considerable shadowing for valence quarks, stronger than for sea 
quarks. We calculate nuclear suppression of DY dilepton production for AuAu collisions as 
function oi xp and impact parameter h and find considerably stronger suppression in AuAu 
collisions than in pAu collisions. Even at large impact parameter, the nucleus-nucleus DY 
cross section is reduced compared to pp as a result of flavor effects. 

We leave for further study the following problems: (i) Since the approximation of a long 
coherence length employed for the DY process at RHIC is satisfied only at > 0.5, one 
must use the LC Green function technique in order to cover the entire range of xp- It can 
also provide a proper interpretation for the Fermilab data |]T|. (ii) Development of the LC 
dipole approach for the DY process in heavy ion collisions is still a challenge; (iii) Nuclear 
modification of the transverse momentum distribution and polarization effects should be also 
calculated for nucleus-nucleus collisions. 
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Appendix A Calculation of gluon shadowing 

Gluon shadowing is given by the shadowing for longitudinal photons as 

^ ^ - AG.{x,Q^) = ' - Aar{x,Q^)- ^^''^ 
and is calculated according to the formulas derived in |2^|. In this appendix, we give some 



details of our calculation. Our starting point is Eq. (90) of [p9f| : 



/poo poc 
d% / dZi / dZ2Q{z2 - Zi)pA{b, Z,)pAib, Z2)T{X, Q\ Z2 - Zi) (A.2) 
J —oo J —oo 



where 

T{x,Q'\Az) = 3? / d\n{aG 
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X [(l-2C-CV^ + C'(3 + C)Ei(C)] 



t sinh(fiAz) 



w 



+ 



In 1 - 



+ 



2t^ t sinh(fiAz) 



+ 



+ 



is given by Eq. (112) of p9| and 

Az = 

n = 
B = 

V 

c 

t 

u 
w 

T2 



Z2 — Zi, 

iB 



"0(1 - «g)z^' 



6^ - iaG(l - OLG)yCeffpA: 



ixitinAz, 
B 

tcosh(fiAz) + sinh(fiAz), 

(1 + t^) sinh(fiA2) + 2t cosh(fiA2) 

(0.65 GeV)^ + acQ"^. 



(A.3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.IO) 
(A.ll) 
(A.12) 



The limits for the ac-integration are x < ac < 0.1, with ac being the momentum fraction 

ar 

«.(Q') 



of the gluon relative to its parent quark. We use a running coupling constant ||59 

An 



q 1 / Q!+a25GeV^ 
(200MeV)2 



(A.13) 



with freezing at low scales. In ( [A.3|) , the gluon-gluon-nucleon cross section is parameterized 
in the form 

a^^{p,x)=C,ff{x)p\ (A.14) 

Note that (Tqq is sampled at the energy x = x/ac- For ac ^ x, x can become greater than 
0.1 and the dipole formulation is no longer valid. To overcome this problem, we employ the 
prescription 

a; = min(x/aG, 0.1). (A. 15) 

The parameter Ce// is then determined from the asymptotic condition 
/ d%d^p |%g(p)|' (1 - exp{-lCeff{x)p^TA{b))) 



Jd'p\^,G{p)fC,ff{x)p' 
Jd'bd'pl^Mf (1 -exp(-K(p,^)r^(6))) 



/rfVI^.G(p)r!<(p 



X) 



(A.16) 



38 



80 



60 \ 



40 



20 





T=2.; 






T=5.3 X 10" 


fm 









0.6 



0.5 



T=2.2 fm"; 
T=5.3 X 10 fm"' 




10"^ 



10 



10"^ 



10"^ 



10"' 



10"^ 



10 



10"^ 



10"' 



10"' 



Figure 17: The x-dependence of Ceff for gold at = 20 GeV^. It is 
also shown how much Ceff is suppressed compared to C{p = 0). The 
two curves in each plot correspond to impact parameter b = (solid) and 
very large impact parameter (dashed), respectively. 



where a^g{p,x) is the dipole cross section in the saturation model of ||T^. The LC wave 
function for radiation of a quark from a gluon, including the nonperturbative interaction 
introduced in |2^, reads 

= \, (A.17) 

The choice of C^ff differs from the one made in |29|, where C(p = 0) = dijQQ{p) / dp'^\p=Q 
was employed as the effective C. The prescription ( A. 16 ) is more realistic, because the Cg// 
is determined by those values of p which are most important for shadowing. Since the 
dipole cross section levels off at large separations, C^ff will be lower than C{p = 0). This 
is illustrated in fig. |T^. The a1 ^ in the denominator of ( |A.2D is calculated with Ceff{x) 
(instead of x). 

With a constant nuclear density p^, one can integrate ( |A.2| ) twice by parts. 



7r 
12 



2Rt 



dL {L^ - 12R\L + V{x, Q^ L), 



(A.18) 



with L = 2^R\ - b^. 

The gg-nucleus cross section in the long coherence time limit is calculated from the 
formula 



atJp,x) = 2 d%{l-{l 



a^g{p,x)TA{b)RG{x,Q\b) 
2A 



(A.19) 



where gluon shadowing as function of impact parameter b is given by 

Aa2*^(x,Q2,6) 



RG{x,Q\b) = l- 



TA{b)anx,Q^y 



(A.20) 
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Figure 18: This figure illustrates 
the influence of gluon shadowing 
on the qq-nucleus cross section. 
It shows the qq-nucleus cross sec- 
tion for gold calculated with gluon 
shadowing, Eg. l \A.19i ), divided by 
the corresponding quantity without 
gluon shadowing, i.e. Rg — > 1- 
The solid curve is for x = 10~^ 
(RHIC), while the dashed curve is 
for X = 10-6 (LHC). 



Aaf^(x, Q\ b)=p\ I dz {L - z) T{x, Q^ z), 



Furthermore, in ( |A.19|) Rq is evaluated at the scale 

1 



+ 4GeV^ 



(A.21) 



(A.22) 



Appendix B Calculation of the DY transverse momen- 
tum distribution 

The differential DY cross section is expressed as a four-fold Fourier integral (|T 
da{qp 7*X) 1 



d In ad^qx 



(2vr; 



2 , d^pid^p2 exp[igT ■ (pi - P2)]^*.9(a, Pi)^7*q(a, P2) 



^ o {^OT-(«Pl) + ^W-(«P2) - CTqqWl - P2))} 



(B.l) 



where 



Up 



27^2 1 Ko(r/pi)Ko(r/p2) 



P1P2 



*:i(a,pi)^^;*o(«,p2 



a 



em ,^2 



■ 7'q 



M^(l-a)^Ko(r/pi)Ko(r/p2) . 



(B.2) 
(B.3) 



and rf = (1 — a)M +a m^. The Fourier-integral is inconvenient for numerical calculations, 
but one can perform three of the integrations analytically for arbitrary cr^(ap). 
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Consider the Kg-part first. With help of the relation 



Ko(r/p) = ^ 



d% 



(B.4) 



one finds 
da{qp 



1*X) 



dhiad?qT 



Off. 



Kq— part 



27r2 



[m)a^ + (1 - af] j d'^pid'^p^ 



d% dH2 
27r 27r 



gi9T (pi-p2) g-i'i-pi gi'2-P2 

2 + ^i5-(«P2) - o^^5-(tt(pi - P2))} • 



(B.5) 



Note that the term in the curly brackets consists of three contributions, which depend either 
only on pi or on p2 or on the difference pi — P2- Thus, the integral (p.5| ) can be split into 
three terms. In the integral that arises from the (^^(api) part, the p2-integration is trivially 

performed and leads to a two dimensional delta-function S^^^qr — h) ■ This makes it possible 
to perform also the I2 integration. The integration over li gives just the MacDonald function 
Kq (p3.4|) . Thus one is left with a two-fold integration over pi. Provided the dipole cross 
section depends only on the modulus of p, one can use the relation 



1 

2^ 



(B.6) 



to perform one more integration. Here, Jq is a Bessel function of first kind. The contribution 
arising from (T^(q;p2) is calculated in exactly the same way. For the (T^(a(pi — p2))-part one 

has to introduce the auxiliary variable d = pi — p2, before the procedure described above 
can be applied. 

The Ki part is calculated in a similar way. Note that 



1 d 

r] dp 



KoM- 



(B.7) 



The Ki-part reads 

da{qp — > 7*X) 



dlnad'^qx 



Ki— part 



[1 



T (27r)2 / 



d'^Pid^P2 



d\ d% 
27r 27r 



X 



X 



(/2 + ^2) (;2 + ^2 



-h ■ h 



a^^{ap2)-a^^{a{p,-p2))]. (B.8) 



Like the Kg-part, the complete integral (|B.8D is split into three pieces, corresponding to the 
three terms in the curly brackets. Again, one integration over p is immediately performed. 
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leading to 5-functions, which allows one to do one integration over /. With the second l- 
integration, one recovers the MacDonald function Ki via ( |B.7| ). For the Ki-part, one also 
needs the relation 

h{z) = -^Jo(^)- (B.9) 

Although the calculation is slightly more cumbersome for the (T^(a(pi — p2))-part, all cal- 
culations are easily performed. 
Finally, one finds 



dlnad'^qT 



27r2 

+ [1 + (!-«)'] 



[mja^ + (1 - af] 



2-^1 ~ l~-^2 



gT^ + r/2 =^ 24' 



with 



X3 



(irrJo(grr)Ko(?7r)cr^(Q;r) 



drr 3o{qTr)Ki{T]r)agg{ar) 



drrJi {qTr)Ki {rir)a Jar) . 



(B.IO) 

(B.ll) 
(B.12) 
(B.13) 



The remaining integrals are evaluated numerically with the Numerical Recipes pDI routines. 



References 

[1] D. M. Aide et al, Phys. Rev. Lett. 64, 2479 (1990). 

[2] M. B. Johnson et al. [FNAL E772 Collaboration], Phys. Rev. Lett. 86, 4483 (2001) 
||hep-ex/0010051|| . 



[3] M. B. Johnson et al., Phys. Rev. C 65, 025203 (2002) ||hep-ph/0105T95|] . 



[4] K. J. Eskola, V. J. Kolhinen and P. V. Ruuskanen, Nucl. Phys. B 535, 351 (1998) 
||hep-ph/9802"350|] ; K. J. Eskola, V. J. Kolhinen and C. A. Salgado, Eur. Phys. J. C 9, 
61 (1999) |Piep-ph/98072g^ . 



[5] L. V. Gribov, E. M. Levin and M. G. Ryskin, Nucl. Phys. B 188, 555 (1981); Phys. 
Rept. 100, 1 (1983). 

[6] P. L. McGaughey, J. M. Moss and J. C. Peng, Ann. Rev. Nucl. Part. Sci. 49, 217 (1999) 
||hep-ph/ 99054091 . 



42 



[7] M. Luo, J. w. Qiu and G. Sterman, Phys. Rev. D 49, 4493 (1994); Phys. Rev. D 
50, 1951 (1994); X. f. Guo, X. f. Zhang and W. Zhu, Phys. Lett. B 476, 316 (2000) 
arXiv:hep-ph/991132l . 



K. J. Eskola, V. J. Kolhinen, P. V. Ruuskanen and R. L. Thews, |hep-ph/0 108093 



B. Z. KopeUovich, proc. of the workshop Hirschegg '95: Dynamical Properties of 
Hadrons in Nuclear Matter, Hirschegg January 16-21, 1995, ed. by H. Feldmeyer and 
W. Ndrenberg, Darmstadt, 1995, p. 102 ( |hep-ph/9609385|) . 

S. D. Drell and T. Yan, Phys. Rev. Lett. 25, 316 (1970) [Erratum-ibid. 25, 902 (1970)]. 
S. J. Brodsky, A. Hebecker and E. Quack, Phys. Rev. D 55, 2584 (1997) [|Eepq 



ph/9609384|| . 



B. Z. Kopeliovich, A. Schafer and A. V. Tarasov, Phys. Rev. C 59, 1609 (1999), extended 
version in Piep-ph/980837^ . 

B. Z. Kopeliovich, J. Raufeisen and A. V. Tarasov, Phys. Lett. B 503, 91 (2001) [|hep^ 
ph/0012035||; J. Raufeisen, J. C. Peng and G. C. Nayak, Phys. Rev. D 66, 034024 (2002) 



|arXiv:hep-ph/0204095|] . 



B. Z. Kopeliovich, J. Raufeisen and A. V. Tarasov, "The color dipole approach to the 
Drell- Yan process in pA collisions," [hep-ph/0104155| . 



S. J. Brodsky, H. C. Pauli and S. S. Pinsky, Phys. Rept. 301, 299 (1998) [ |arXiv:hep- 
ph/9705477|| . 



A. B. Zamolodchikov, B. Z. Kopeliovich and L. I. Lapidus, JETP Lett. 33, 595 (1981) 
[Pisma Zh. Eksp. Teor. Fiz. 33, 612 (1981)]. 

K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59, 014017 (1999) ||hep-ph/98075T3 
Phys. Rev. D 60, 114023 (1999) ||hep-ph/ 990335^ 



M. B. Johnson, B. Z. Kopeliovich and A. V. Tarasov, Phys. Rev. C 63, 035203 (2001) 
||hep-ph/ 00063^ ]. 



G. T. Bodwin, S. J. Brodsky and G. P. Lepage, Phys. Rev. D 39, 3287 (1989). 

G. Bertsch, S. J. Brodsky, A. S. Goldhaber and J. F. Gunion, Phys. Rev. Lett. 47, 297 
(1981). 

C. S. Lam and W.-K. Tung, Phys. Rev. D 18, 2447 (1978); Phys. Rev. D 21, 2712 
(1980). 

R. J. Glauber "High Energy Collision theory" in Lectures in Theoretical Physics, vol. 1, 
W. E. Brittin and L. G. Duham (eds.) Interscience, New York, 1959. 



43 



[23 



[24 

[25' 
[26 



[27 



[28 



[29 
[30 



[31 



[32 
[33 
[34 
[35; 
[36 

[37; 

[3 



Y. V. Kovchegov and A. H. Mueller, Nucl. Phys. B 529, 451 (1998) [|arXiv:liep- 
ph/980244q| . 



S. J. Brodsky, P. Hoyer, N. Marchal, S. Peigne and F. Sannino, Phys. Rev. D 65, 114025 
(2002) [|arXiv:liep-pli/0104291|| . 



U. A. Wiedemann, Nucl. Phys. B 582, 409 (2000) [ |arXiv:hep-ph/ 000302 1|| . 
A. H. Mueller, Nucl. Phys. B 335, 115 (1990); Nucl. Phys. B 558, 285 (1999) [|Ei^ 



ph/9904404|| . 



O. V. Kancheli, Print-75-0378 (TBILISI) In *Caneschi, L. (ed.): Regge theory of low- 
p(T) hadronic interactions* 312-315. (JETP Lett. 18 (1973) 274-277). (Pisma Zh. Eksp. 
Teor. Fiz. 18 (1973) 465-469). (see Book Index). 

B. Blattel, G. Baym, L. L. Frankfurt and M. Strikman, Phys. Rev. Lett. 70, 896 (1993); 
L. Frankfurt, A. Radyushkin and M. Strikman, Phys. Rev. D 55, 98 (1997) ||arXiv:hep^ 
ph/96 102741] . 



B. Z. Kopeliovich, A. Schafer and A. V. Tarasov, Phys. Rev. D 62, 054022 (2000) 
||hep-ph/ 99082151 ]. 



A. L. Ayala, M. B. Gay Ducati and E. M. Levin, Nucl. Phys. B 493, 305 (1997) [|h^ 
ph/9604383|| . 



B. Kopeliovich, A. Tarasov and J. Hiifner, Nucl. Phys. A 696, 669 (2001) [ |hep- 
ph/0104256^ . 



B. Z. Kopeliovich, J. Nemchik, A. Schafer and A. V. Tarasov, Phys. Rev. C 65, 035201 
(2002) ||hep-ph/ 01072271 . 

V. N. Gribov, Sov. Phys. JETP 29, 483 (1969) [Zh. Eksp. Teor. Fiz. 56, 892 (1969)]; 
Sov. Phys. JETP 30, 709 (1970) [Zh. Eksp. Teor. Fiz. 57, 1306 (1970)]. 

I. Balitsky, Nucl. Phys. B 463, 99 (1996) [ |arXiv:hep-ph/9509348|| . Y. V. Kovchegov, 
Phys. Rev. D 60, 034008 (1999) ||arXiv:hep-ph/9901281|| . 



C. W. De Jager, H. De Vries and C. De Vries, Atom. Data Nucl. Data Tabl. 36, 495 
(1987). 

T. Schafer and E. V. Shuryak, Rev. Mod. Phys. 70, 323 (1998) | |hep-ph/9610"45T| . 

V. M. Braun, P. Gornicki, L. Mankiewicz and A. Schafer, Phys. Lett. B 302, 291 (1993). 

M. D'Elia, A. Di Giacomo and E. Meggiolaro, Phys. Lett. B 408, 315 (1997) [|hep^ 



I lat/9705032i . 

[39] B. Z. Kopeliovich, J. Raufeisen and A. V. Tarasov, Phys. Rev. C 62, 035204 (2000) 
||hep-ph/000313q . 



44 



H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12, 375 (2000) ||liep-pli/990328^ . 

H. Plothow-Bescli, "The Parton distribution function library," Int. J. Mod. Phys. A 10, 
2901 (1995); "PDFLIB: Proton, Pion and Photon Parton Density Functions, Parton 
Density Functions of the Nucleus and Calculations", User's Manual - Version 8.04, 
W5051 PDFLIB, 2000.04.17, CERN-PPE. 

J. Raufeisen, Ph.D. thesis "QCD coherence effects in high energy reactions with nuclei," 
|hep-ph/0009"358| . 

A. Drees, Nucl. Phys. A 698, 331 (2002) [|nucl-ex/0 1050191] . 

B. Z. Kopeliovich and F. Niedermayer, Nuclear screening in J/\I/ and Drell-Yan pair 
production, JINR-E2-84-834, Dubna 1984, a scanned version in KEK library: 

|http: / / www-lib.kek. jp/ cgi-bin/img_index?85041 13 . 

J.-C. Peng, private communication. 

D. M. Aide et al, Phys. Rev. Lett. 66, 2285 (1991). 

D. M. Kaplan et al, Phys. Rev. Lett. 40, 435 (1978). 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994) ||hep-ph/9309289|] ; 
Phys. Rev. D 49, 3352 (1994) ||hep-ph/9311205|] ; Phys. Rev. D 50, 2225 (1994) 
ph/9402335| |. 



J. Dolejsi, J. Hiifner and B. Z. Kopeliovich, Phys. Lett. B 312, 235 (1993) [fiep^ 
ph/9305238|] . 



E. Mirkes, Nucl. Phys. B 387, 3 (1992); A. Brandenburg, O. Nachtmann and E. Mirkes, 
Z. Phys. C 60, 697 (1993). 

S. Falciano et al. [NAlO Collaboration], Z. Phys. C 31, 513 (1986); J. S. Conway et al, 
Phys. Rev. D 39, 92 (1989). 



R. J. Fries, A. Schafer, E. Stein and B. Miiller, Nucl. Phys. B 582, 537 (2000) [^ep 



ph/0002074 ; R. J. Fries, B. Miiller, A. Schafer and E. Stein, Phys. Rev. Lett. 83, 4261 



(1999) [|iep-ph/9907567|| . 

J. C. Collins, D. E. Soper and G. Sterman, Nucl. Phys. B 250, 199 (1985). 

M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 410, 337 (1997). 

J. C. Peng et al. [E866/NuSea Collaboration], Phys. Rev. D 58, 092004 (1998) 



ph/9804288| . 



[56] B. Z. Kopeliovich, J. Raufeisen and A. V. Tarasov, Phys. Lett. B 440, 151 (1998) [fiep^ 



I ph/980721l| ; J. Raufeisen, A. V. Tarasov and O. O. Voskresenskaya, Eur. Phys. J. A 
5, 173 (1999) ||hep-ph/9812398|] . 



45 



[57] M. Arneodo, Phys. Kept. 240, 301 (1994). 

[58] M. Gliick, E. Reya and A. Vogt, Eur. Phys. J. C 5, 461 (1998) ||hep-ph/9806l0^] . 

[59] Y. L. Dokshitzer, G. Marchesini and B. R. Webber, Nucl. Phys. B 469, 93 (1996) 
||hep-ph/951233q] . 

[60] W. H. Press, S. A. Teukolsky, W. T. Vetterhng and B. P. Flannery, "Numerical Recipes 
in C" , 2^^*^ edition, Cambridge University Press, New York, 1995. 



46 



